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upper  ocean.  Recent  developments  in  several  technologies  and  modeling  capabilities  have 
enabled  the  investigation  of  a  variety  of  fundamental  and  applied  problems  related  to 
upper  ocean  physics,  chemistry,  and  light  propagation  and  utilization  in  the  dynamic 
near-surface  ocean.  The  purpose  here  is  to  provide  background  for  and  an  introduction 
to  a  collection  of  papers  devoted  to  new  technologies  and  observational  results  as  well 
as  model  simulations,  which  are  facilitating  new  insights  into  optical  variability  and 
light  propagation  in  the  ocean  as  they  are  affected  by  changing  atmospheric  and 
oceanic  conditions. 

Citation:  Dickey,  T.,  et  al.  (2012),  Introduction  to  special  section  on  Recent  Advances  in  the  Study  of  Optical  Variability 
in  the  Near-Surface  and  Upper  Ocean,  J.  Geophys.  Res.,  117,  C00H20,  doi:l 0.1 029/201 2 JC007964. 


1.  Introduction 

[2]  Many  of  the  papers  in  this  special  section  were  stimu¬ 
lated  by  research  conducted  during  the  Radiance  in  a 
Dynamic  Ocean  (RaDyO)  program  | Dickey  et  al,  2011] 
while  others  are  from  related  studies.  One  of  the  primary 
goals  of  RaDyO  was  to  study  the  problem  of  underwater 
visibility  and  imaging  of  objects  across  the  air-sea  interface 
and  of  objects  residing  within  the  upper  ocean  as  affected  by 
a  broad  range  of  environmental  conditions  [e.g.,  Duntley, 
1963;  Mullamaa,  1975;  Preisendorfer,  1976;  Davies- 
Collev,  1988;  Schippnick,  1988;  McLean  and  Freeman, 
1996;  Jaffe  et  al.,  2001;  Zaneveld  and  Pegau,  2003;  Dolin 
et  al.,  2006;  Doron  et  al.,  2007;  Dolin  and  Luchinin,  2008; 
Levin  et  al.,  2008].  This  specific  problem  and  other  related 
problems,  as  described  below,  involve  meteorological, 
physical,  optical,  biological,  and  chemical  processes  acting 
over  time  and  space  scales  as  short  as  milliseconds  and  mil¬ 
limeters.  For  example,  atmospheric  conditions  involving 
clouds,  fog,  aerosols,  and  sun  angle  affect  incident  solar 
radiation  impinging  on  the  ocean’s  surface.  The  initial  fate  of 
the  impinging  light  is  controlled  by  the  air-sea  interface’s 
topography,  as  modulated  by  surface  capillary  and  gravity 
waves,  and  internal  gravity  waves  [e.g.,  Schenck,  1957; 
Snyder  and  Dera,  1970;  Dera  and  Olszewski,  1978;  Stramski 
and  Dera,  1988;  Dera,  1992;  Dera  and  Stramski,  1993; 
McLean  and  Freeman,  1996;  Stramska  and  Dickey,  1998; 
Zaneveld  et  al.,  2001a,  2001b;  Wijesekera  et  al.,  2005;  You 
et  al.,  2009,  2010],  and  as  affected  by  the  presence  of 
surfactants  [e.g.,  Wurl  et  al.,  2009]  and  bubbles  [e.g., 
Stramski,  1994;  Zhang  et  al.,  1998;  Stramski  and  Tegowski, 
2001;  Terrill  et  al.,  2001;  Wolk  et  al.,  2002;  Terrill  and 
Lewis,  2004;  Zhang  et  al.,  2011],  Light  propagation  within 
and  exiting  the  upper  ocean  is  also  affected  by  organisms, 
particulate  and  dissolved  materials  including  colored  dis¬ 
solved  organic  matter,  physical,  chemical, and  biological 
processes,  inherent  optical  properties  (lOPs),  and  apparent 
optical  properties  (AOPs)  along  with  characteristics  of  the 
molecular  and  turbulent  boundary  layers  on  both  sides  of  the 
air-sea  interface.  These  layers  are  dynamic  and  affected  by 
wind-forcing,  waves,  stratification,  chemical  and  particulate 
composition  of  seawater,  optical  properties,  bubbles,  foam, 
and  whitecaps  (W.  K.  Melville  et  al.,  The  impact  of  surface 
wave  breaking  on  imaging  through  the  ocean  surface,  man¬ 
uscript  in  preparation,  2012). 

[3]  There  are  several  other  relevant  problems  that  relate  to 
this  research  and  require  similar  data  sets  and  modeling 
approaches.  These  include  the  fundamental  problem  of 
radiative  transfer  (e.g.,  light  propagation  across  the  air-sea 


interface  and  within  the  upper  layer  [e.g.,  Kattawar  et  al., 
1976;  Shifrin,  1988;  Kattawar  and  Adams,  1989;  Dera, 
1992;  Kirk,  1994;  Mobley,  1994;  Walker,  1994;  Thomas 
and  Stamnes,  1999;  Zhai  et  al.,  2008a,  2008b]);  marine  ani¬ 
mal  vision  [e.g.,  McFarland  andLoew,  1983],  phytoplankton 
responses  to  fluctuating  light  fields  [e.g.,  Frechette  and 
Legendre,  1978;  Walsh  and  Legendre,  1983;  Legendre  et  al, 
1986;  Queguiner  and  Legendre,  1986;  Stramski  et  al.,  1993; 
Kirk,  1994],  air-sea  exchanges  of  momentum,  heat,  and  gases 
[e.g.,  Thorpe,  1982;  Wanninkhof  1992;  Woolf,  1997,  2005; 
Weller  et  al.,  1998;  Wanninkhof  and  McGillis,  1999; 
Wanninkhof  et  al.,  2004;  Yu  and  Weller,  2007;  Berry  and 
Kent,  2009;  Brooks  et  al,  2009],  and  the  interpretation  of 
radiometric  data  collected  in  the  ocean  from  various  plat¬ 
forms  including  imagers  mounted  on  aircraft  and  satellites 
[e.g.,  Fraser  et  al,  1980;  McLean  and  Freeman,  1996; 
Zaneveld  et  al,  2001a,  2001b;  Zibordi  et  al.,  2004; 
Twardowski  et  al.,  2007a;  Zaneveld  et  al,  2005;  Lewis, 
2008;  Gernez  and  Antoine,  2009], 

[4]  Fluctuations  in  the  underwater  light  field  caused  by 
surface  waves  have  been  a  subject  of  earlier  theoretical  and 
experimental  research.  The  models  of  Schenck  [1957]  and 
Snyder  and  Dera  [1970]  identified  mechanisms  for  produc¬ 
ing  flashing  light,  which  are  associated  with  fluctuations  in 
sea  surface  topography,  namely  its  elevation,  surface  slope, 
and  surface  curvature.  These  physical  processes  are  highly 
complex  and  can  invalidate  simplifying  assumptions;  thus, 
the  solution  of  the  theoretical  problem  is  challenging.  Light 
fluctuations  just  beneath  the  ocean  surface  are  most  intense 
under  sunny  skies  accompanied  by  relatively  weak  winds 
and  low  sea  states  when  wave  focusing  of  direct  solar  rays 
dominates.  The  intensity  and  high-frequency  content  of 
these  fluctuations  decrease  rapidly  with  depth  within  the  top 
several  meters  of  the  ocean  [e.g.,  Dera  and  Gordon,  1968; 
Snyder  and  Dera,  1970;  Fraser  et  al.,  1980;  Dera  and 
Stramski,  1986]. 

[5]  Specially  designed  radiometers  and  sampling  strate¬ 
gies  are  required  to  fully  resolve  the  maximum  focusing  at 
near-surface  depths  (down  to  about  5-10  m  in  clear  waters), 
because  these  effects  are  often  characterized  by  the  presence 
of  light  pulses  of  very  high  amplitude  and  short  duration 
[Dera  and  Olszewski,  1978;  Dera  and  Stramski,  1986; 
Gernez  et  al.,  2011;  Darecki  et  al,  2011].  Most  experi¬ 
mental  work  on  underwater  light  fluctuations  to  date  has 
been  best  suited  to  examine  weaker  fluctuations  of  light  that 
occur  at  depths  below  the  top  few  meters  of  the  ocean.  At 
these  depths,  even  under  sunny  conditions,  the  fluctuations 
in  instantaneous  irradiance  generally  do  not  exceed  the  mean 
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irradiance  by  more  than  50%,  and  are  often  within  10-20% 
of  the  mean  irradiance  level.  The  book  by  Walker  [1994] 
reviews  early  relevant  studies,  including  several  from  the 
former  Soviet  Union.  Examples  of  more  recent  work  include 
Stramska  and  Dickey  [1998]  and  Gernez  and  Antoine 
[2009],  who  noted  that  if  skies  are  overcast  and  surface 
illumination  is  diffuse,  underwater  light  fluctuations  are 
weak  throughout  the  entire  water  column  [see  also  Dera  and 
Olszewski ,  1967;  Stramski  et  al.,  1992],  These  fluctuations 
consist  primarily  of  low-frequency  components  due  to  sur¬ 
face  displacements  associated  with  significant  surface  waves 
and/or  swell. 

[6]  Interacting  processes  necessitate  novel  methods  and 
models  for  imaging  and  other  applications.  A  major  obstacle 
to  this  line  of  research  has  been  the  inability  of  sampling  and 
computational  methods  and  hardware  to  sufficiently  resolve 
high  temporal  and  spatial  scale  variabilities,  which  can  have 
very  large  dynamic  ranges  (up  to  at  least  6  to  8  orders  of 
magnitude)  for  such  a  diverse  set  of  physical  and  optical 
variables.  Several  previous  studies  have  used  instrumenta¬ 
tion  designed  specifically  to  measure  the  fast  fluctuations 
and  associated  strongest  wave  focusing  effects  occurring  at 
near-surface  depths  under  sunny  conditions  [Dera  and 
Olszewski,  1978;  Dera  and  Stramski,  1986;  Stramski, 
1986;  Dera  and  Stramski,  1993].  However,  this  work  has 
been  mostly  limited  to  a  single  light  wavelength,  typically 
selected  for  maximum  transparency  in  water.  Focusing  has 
been  found  to  be  most  pronounced  under  weak  winds  (2- 
5  m/s)  and  thus  slightly  undulating  sea  surfaces,  clear 
atmospheres  with  low  surface  irradiance  diffuseness  of  less 
than  40%,  and  relatively  high  sun  elevation  with  solar  zenith 
angle  generally  less  than  40°  [Dera  and  Stramski,  1986; 
Stramski,  1986].  Under  such  conditions,  high-amplitude 
pulses  of  focused  light,  exceeding  the  mean  level  of  irradi¬ 
ance  by  more  than  50%,  have  been  shown  to  occur  at  fre¬ 
quencies  higher  than  100  per  minute  at  near-surface  depths. 
These  strong  focusing  events,  referred  to  as  light  flashes, 
have  short  durations  on  the  order  of  10  ms.  Reviews  of 
general  features  of  wave-induced  light  fluctuations  are  pro¬ 
vided  in  papers  by  Stramski  et  al.  [1992],  Gernez  and 
Antoine  [2009],  Gernez  et  al.  [2011],  and  Darecki  et  al. 
[2011], 

[7]  Previous  work  on  light  fluctuations  has  been  unable  to 
include  some  of  the  more  important  physical  variables  in 
sufficient  detail,  especially  wave  characteristics.  RaDyO 
researchers  are  beginning  to  report  significant  advances  in 
terms  of  linking  light  fluctuations  to  surface  waves,  because 
specially  designed  instrumentation  and  experiments  have 
been  utilized  recently  to  obtain  relevant  data  sets  for  devel¬ 
oping  and  parameterizing  more  realistic  models  that  couple 
radiative  transfer,  surface  waves,  and  other  physical 
dynamics. 

[8]  It  should  be  noted  that  there  has  been  an  increased 
awareness  of  the  multiplicity  of  interdisciplinary  processes 
that  must  be  considered  in  studying  problems  such  as  optical 
variability  of  the  ocean.  For  example,  bubbles,  which  are 
ubiquitous  in  the  world  ocean,  play  important  roles  for  a  host 
of  processes.  The  problem  of  estimating  the  number  of 
bubbles  and  their  size  distributions  is  daunting.  MacIntyre 
[1972]  forwarded  an  estimate  that  some  10 18  to  1020  white 
cap  bubbles  break  (and  fonn)  per  second  over  the  world 
ocean.  More  recently,  Bird  et  al.  [2010]  have  suggested  that 


the  number  of  bubbles  may  be  even  greater  based  on  their 
work  on  daughter  bubble  cascades  (creations  of  toroidal 
rings  of  small  bubbles  upon  rupture).  The  importance  of 
bubbles  for  optics  of  the  ocean  has  been  recognized  by  a  few 
investigators  [e.g.,  Stramski,  1994;  Zhang  et  al.,  1998; 
Stramski  and  Tegowski,  2001;  Terrill  et  al.,  2001;  Terrill 
and  Lewis,  2004;  Zhang  et  al.,  2011].  In  fact,  it  has  been 
suggested  that  a  large  portion  of  upwelled  radiance  from  the 
ocean  surface  may  be  attributed  to  bubbles  and  that  ocean 
color  shifts  (toward  green  wavelengths)  can  affect  ocean 
remote  sensing  algorithms  for  chlorophyll  a  [e.g.,  Zhang 
et  al.,  1998].  Bubble  size  distributions  and  organic  coat¬ 
ings  of  bubbles  play  important  roles  as  well  [e.g.,  Zhang 
et  al.,  1998;  Stramski  and  Tegowski,  2001;  Zhang  et  al., 
2011].  The  recognition  of  the  potential  optical  effects  of 
bubbles  has  led  to  the  development  of  new  optical  and 
acoustical  instrumentation  with  the  expressed  purposes  of 
quantifying  bubble  effects  on  light  in  the  ocean;  some  of 
these  are  described  in  this  special  section  [Czerski  et  al., 
2011;  Vagle  et  al.,  2012;  Twardowski  et  al.,  2012].  Impor¬ 
tantly,  overall  sampling  of  ocean  optical  variability  has  been 
advanced  by  the  utilization  of  multiple  platforms  [e.g., 
Dickey  and  Bidigare,  2005;  Dickey  et  al.,  2006],  which 
facilitate  the  concurrent  collection  of  a  plethora  of  interdis¬ 
ciplinary  data  sets  on  unprecedented  time  and  space  scales. 
The  RaDyO  program  was  also  devoted  to  modeling  optical 
variability  of  the  surface  boundary  layer  and  upper  ocean, 
likewise  on  unprecedented  time  and  space  scales.  Some  of 
this  modeling  work  is  introduced  in  section  5. 

[9]  This  paper  concludes  by  briefly  introducing  contribu¬ 
tions  to  this  special  section  categorically  as  follows:  (1) 
observational  papers  based  primarily  on  results  from  the 
RaDyO  program,  (2)  modeling  papers  done  in  conjunction 
with  RaDyO,  and  (3)  papers  relevant  to  the  theme  of  the 
special  section  that  were  submitted  by  authors  independently 
of  RaDyO.  Finally,  interested  readers  are  directed  to  a  recent 
survey  article  by  Dickey  et  al.  [2011],  which  briefly  provides 
historical,  introductory,  and  summary  information  relevant 
to  this  special  section. 

2.  Overview  of  the  RaDyO  Observational 
Program 

2.1.  Sites  of  RaDyO  Observations 

[10]  Fulfillment  of  the  objectives  of  RaDyO  required  that  a 
wide  range  of  environmental  conditions  be  observed  for  the 
development  of  models  that  would  be  useful  for  diverse 
optical,  meteorological,  and  physical  situations  [Dickey 
et  al.,  2011].  The  first  RaDyO  experiment  was  conducted 
using  the  Scripps  Institution  of  Oceanography  (SIO)  pier  6- 
28  January  2008.  This  work  was  primarily  intended  to  test 
new  instrumentation  in  a  setting  with  good  access  to  shore- 
based  technical  support  and  to  develop  integration  of  the 
diverse  methodologies.  However,  some  scientific  results 
have  been  produced  from  this  experiment  as  well.  The  sec¬ 
ond  experiment  was  conducted  in  the  Santa  Barbara  Channel 
(SBC)  during  the  period  of  September  3-25,  2008.  This 
setting  was  selected  because  it  was  expected  to  provide  a 
relatively  benign  wind-wave  regime  and  easy  access  to  shore 
since  several  new  or  prototype  instrumentation  systems  were 
being  utilized.  The  third  and  final  RaDyO  experiment  took 
place  south  of  the  Big  Island  of  Hawaii  from  24  August 
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Figure  1.  Photographs  of  the  R/V  Kilo  Mo  ana.  Lil  KM,  R/P  FLIP,  Odyssey’  AUV,  and  REMUS  AUV. 


through  September  15,  2009.  This  location  was  selected 
because  of  its  climatologically  high,  persistent  wind  speeds 
and  relatively  large  sea  states,  its  optically  clear  waters,  and 
its  open  ocean  character  with  good  access  to  the  deep-water 
port  of  Honolulu.  The  specific  area  of  operations,  which  lies 
in  the  North  Equatorial  Current  (NEC),  is  characterized  by 
relatively  modest  mesoscale  eddy  activity  and  minor  island 
effects  as  compared  with  areas  directly  in  the  wind  wake  of 
the  Hawaiian  archipelago  [e.g.,  see  Dickey  et  al.,  2008].  The 
field  experiments  are  described  in  more  detail  below  and  on 
the  Website  http://www.opl.ucsb.edu/radyo/. 

2.2.  Platforms 

[1 1]  The  multifaceted,  multiscale  nature  of  RaDyO  pro¬ 
blems  necessitated  a  multiple  platform,  multidisciplinary 
sampling  approach  as  has  been  discussed  in  detail  in  earlier 
references.  The  primary  platforms  used  for  the  RaDyO  Santa 
Barbara  Channel  experiment  included:  R/P  FLIP,  R/V  Kilo 
Moana,  an  Odyssey  autonomous  underwater  vehicle  (AUV) 
[Moline  et  al.,  2012]  (Figure  1),  a  REMUS  AUV  [Moline 
et  al.,  2005,  2007,  2010,  2012;  Robbins  et  al.,  2006],  and 
a  surfactant  skimmer  called  the  Lil  KM  (Figure  1).  Also,  a 
small  aircraft  equipped  with  lidar  instrumentation  made 
measurements  over  the  site  [Reineman  et  al,  2009].  Because 
of  anticipated  greater  sea  states,  the  Odyssey  AUV  and  the 
Lil  KM  were  not  utilized  and  the  REMUS  AUV  was 
deployed  on  only  a  limited  basis  for  the  Hawaiian  experi¬ 
ment.  In  addition,  satellite-deployed  sensors  were  utilized  to 
obtain  wind,  sea  surface  temperature,  color,  and  surface 
roughness  (SAR  for  Santa  Barbara  experiment  only)  data. 
The  platforms  used  for  the  field  experiments  are  briefly 
described  below. 

2.2.1.  RIP  FLIP 

[12]  R/P  FLIP  (Floating  Instrument  Platform,  Figure  2), 
originally  launched  on  2  June  1962,  was  designed  to  be  a 
highly  stable  platform  for  oceanographic  research  in  the 
areas  of  acoustics,  seismic  waves,  ocean  waves,  and  marine 
biology  [Fisher  and  Spiess,  1963].  However,  it  has  since 


been  used  for  a  variety  of  meteorological,  air-sea  interaction, 
and  physical  oceanographic  experiments  [e.g.,  Siegel  and 
Dickey,  1986;  Smith  et  al.,  1987;  Pinkel  et  al.,  1991; 
Plueddemann  et  al.,  1996;  Smith,  1998;  Friehe  et  al.,  2001; 
Rainville  and  Pinkel,  2006]  and  for  optical  measurements  in 
1982  [Siegel  and  Dickey,  1987;  Dickey,  1991],  R/P  FLIP  is 
essentially  a  large  manned  spar  buoy  when  oriented  verti¬ 
cally  (note  that  it  is  towed  horizontally  to  experiment  sites). 
R/P  FLIP  measures  approximately  108  m  in  length  with 
~90  m  of  its  hull  extending  below  the  water  line  when  on 
station.  The  diameter  of  R/P  FLIP’S  hull  is  6.5  m  from  the  49 
to  91  m  depth;  the  hull  tapers  to  4  m  at  the  20  m  depth  (see 
schematic  in  Figure  3  of  Fisher  and  Spiess  [1963]).  This 
design  allows  R/P  FLIP,  in  free  drift  mode,  to  be  minimally 
responsive  to  wave  motion  (<10%  of  surface  motion,  reso¬ 
nant  period  is  ~25  s  and  thus  displaced  away  from  the 
spectral  window  of  greatest  wave  motion  (about  5-18  s))  as 
documented  by  Fisher  and  Spiess  [1963],  Rudnick  [1967], 
Rudnick  and  Hasse  [1971],  and  Smith  and  Rieder  [1997]. 
For  example,  in  ten  meter  waves  R/P  FLIP’S  total  vertical 
motion  has  been  recorded  to  be  less  than  1  m  (see  http:// 
www.mpl.ucsd.edu/resources/flip.intro.html  for  a  bibliogra¬ 
phy  and  details  concerning  R/P  FLIP). 

[13]  For  the  Santa  Barbara  Channel  RaDyO  experiment, 
R/P  FLIP  (Figure  2a)  was  moored  in  place  at  34.2053°  N, 
119.6288°  W  in  water  of  168  m  depth  using  a  3-point 
mooring  [e.g.,  see  Bronson,  1971],  Some  tilting  of  R/P  FLIP 
was  apparent  due  to  currents  and  winds.  These  conditions 
were  apparently  accompanied  by  a  measurable  oscillation  in 
R/P  FLIP’S  heading  and  depth  with  an  approximately  2  min 
period.  Instruments  were  deployed  from  R/P  FLIP’S  booms 
to  minimize  the  influence  of  flow  distortion  and  super¬ 
structure  interference  with  measurement  devices  which  were 
either  profiled  or  deployed  at  or  mounted  at  specific  depths. 
Booms  with  trolleys  were  used  to  deploy  instrumentation. 
Also,  R/P  FLIP  necessarily  discharged  wastewaters  period¬ 
ically  and  the  times  of  these  discharges  were  noted  for  pos¬ 
sible  contamination  of  optical  and  chemical  data  sets.  For 
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(a)  R/P  FLIP  •  Santa  Barbara  Experiment  (2008) 
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(b)  R/P  FLIP  -  Hawaii  Experiment  (2009) 
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Figure  2.  Instrumentation  deployed  on  R/P  FLIP  during  the  (a)  Santa  Barbara  and  (b)  Hawaii  experi¬ 
ments  in  September  2008  and  2009,  respectively,  (c)  R/P  FLIP  starboard  boom  during  the  Hawaii  Exper¬ 
iment  in  September  2009.  The  air-sea  flux  package,  orthogonal  scanning  laser  altimeters,  polarimetric  and 
infrared  cameras  are  shown  in  the  upper  right  insert.  The  acoustic  resonators  and  sonar  float  [Vagle  et  al., 
2012]  are  highlighted  in  the  lower  insert,  (d)  Corresponding  instrumented  R/P  FLIP  port  boom,  showing  a 
suite  of  air-sea  interaction  and  subsurface  electro-optical  sensors  used  to  characterize  the  influence  of 
breaking  at  the  ocean  surface  on  oceanic  radiance  and  imaging  (Melville  et  al.,  manuscript  in  preparation, 
2012):  eddy  covariance  meteorological  system,  stereo  infrared  and  visible  cameras,  scanning  and  single¬ 
point  laser  altimeters,  LED  downward  looking  display  for  imaging  through  the  surface  with  an  underwater 
upward  looking  color  camera  at  various  depths.  This  subsurface  profiler  is  also  equipped  with  a  Nortek 
HR  Profiler  and  ADV  for  current  and  turbulence  measurements. 


reference,  R/V  Kilo  Moana  sampled  on  station  about  2  km 
north  of  R/P  FLIP  during  the  SBC  experiment  except  during 
evening  treks  away  from  the  site  to  discharge  its  tanks. 

[14]  The  Odyssey  and  REMUS  AUVs  conducted  their 
mapping  surveys  between  R/P  FLIP  and  R/V  Kilo  Moana  as 
well  as  in  their  vicinities  in  patterns  described  by  Moline 
et  al.  [2012].  Finally,  it  should  be  noted  that  R/P  FLIP 
(Figure  2b)  was  allowed  to  drift  in  response  to  winds  and 
currents  during  the  Hawaii  RaDyO  experiment  because 


mooring  it  in  deep  water  was  not  feasible  for  our  experiment. 
Furthermore,  tilt  effects  were  minimal  in  this  case.  The  R/V 
Kilo  Moana  again  sampled  as  close  to  R/P  FLIP  as  safely 
allowable  (generally  within  2  km)  and  the  REMUS  AUV 
was  occasionally  used  for  spatial  sampling. 

2.2.2.  R/V  Kilo  Moana 

[is]  The  R/V  Kilo  Moana  is  a  Small  Waterplane  Area 
Twin  Hull  (SWATH)  vessel  of  56.7  m  length  and  with  a 
beam  of  26.8  and  a  draft  of  7. 0-7. 6  m  [see  http://www.soest. 
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Figure  2.  (continued) 


hawaii.edu/UMC/KiloMoana.htm  for  details].  The  R/V  Kilo 
Moana,  which  was  designed  to  minimize  adverse  motions 
via  its  design  and  motion  compensation  system  in  moderate 
to  high  sea  states,  was  used  for  a  variety  of  RaDyO  mea¬ 
surements.  In  sea  state  5  with  significant  wave  heights  of 
2. 5-4.0  m,  R/V  Kilo  Moana  ’s  maximum  pitch  is  9.2  deg, 
maximum  heave  is  8.2  m  and  maximum  roll  is  15.2  deg  (all 
measured  at  the  vessel’s  center  of  gravity).  The  instrumen¬ 
tation  deployed  from  R/P  FLIP  and  the  R/V  Kilo  Moana  is 
introduced  below  and  discussed  in  more  detail  in  other 
papers  in  this  special  section. 

2.3.  Measurement  Systems 

[i6]  Several  standard  and  new  technologies  were  critical 
for  meeting  the  objectives  of  RaDyO.  Some  of  the  novel 
technologies  included:  (1)  optical  imaging  systems  capable 
of  sampling  at  several  Hertz  and  with  spatial  resolutions 
down  to  a  few  mm,  (2)  new  imaging  systems  with  high 
scene  dynamic  ranges  to  image  both  direct  and  diffuse 
radiance  fields,  (3)  polarigraphic  optical  measurements,  (4) 


underwater  radiometric  systems  consisting  of  numerous 
radiance  and  irradiance  sensors  for  fast  sampling  of  wave- 
induced  light  fluctuations  at  1  kHz  frequency,  (5)  optical 
instruments  for  measuring  absorption  and  scattering  (some 
at  multiple  scattering  angles)  at  multiple  wavelengths  in  the 
visible  spectrum,  (6)  meteorological  instruments  for  directly 
measuring  momentum  and  heat  fluxes,  (7)  a  small  surface 
skimming  vessel,  Lil  KM,  for  collection  of  microlayer  sur¬ 
factants  and  physical  and  optical  measurements  in  the  upper 
1.7  m  of  the  water  column,  and  (8)  two  autonomous  under¬ 
water  vehicles  for  sampling  physical  and  optical  variability. 

[n]  For  convenience,  the  general  descriptions  of  RaDyO 
measurements  are  subdivided  first  by  platform  and  then  by 
disciplinary  type.  Most  of  the  systems  were  similar  for  the 
SBC  and  Hawaii  experiments,  however  significant  differ¬ 
ences  are  noted  below  for  completeness.  Readers  interested 
in  more  detailed  technical  information  and  the  names  of 
investigators  who  were  responsible  for  specific  measure¬ 
ments  are  directed  to  the  Website  www.opl.ucsb.edu/radyo/ 
and  other  papers  in  this  special  section. 
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2.3.1.  R/V  Kilo  Moana  Measurements 

[is]  Observations  were  made  from  the  R/V  Kilo  Moana 
using  instrumentation  provided  by  the  University  of 
Hawaii’s  instrumentation  pool  and  by  RaDyO  research 
groups.  During  each  cruise,  the  R/V  Kilo  Moana’ s  underway 
measurement  system  collected  a  variety  of  data  at  1  to  10  s 
intervals  depending  on  the  variable  and  final  processed 
data  were  averaged  on  time  scales  from  1  s  to  5  min. 
The  measurement  systems  deployed  from  the  R/V  Kilo 
Moana  during  RaDyO  are  summarized  below,  on  the 
RaDyO  website,  and  at  http://www.soest.hawaii.edu/UMC/ 
KM/ scienceequipment.htm. 

2.3.1. 1.  Meteorological  Measurements 

[19]  The  R/V  Kilo  Moana’s  meteorological  sensors  were 
located  at  20.7  m  above  the  sea  surface  and  data  were  col¬ 
lected  through  the  ship’s  underway  sampling  system.  Vari¬ 
ables  included:  GPS  position,  meteorological  variables  of 
barometric  pressure  (Vaisala),  wind  speed  and  direction 
(RM  Young),  relative  humidity  (MP101A-C),  precipitation 
(Rotronics),  air  temperature,  shortwave  and  longwave  radi¬ 
ation  (Eppley  PSP  and  PIR),  and  surface  photosynthetically 
available  radiation  (PAR). 

2.3. 1.2.  Physical  Measurements 

[20]  Underway  data  were  collected  through  the  ship’s 
intake,  which  was  located  8  m  below  the  sea  surface,  for 
obtaining  uncontaminated  seawater.  These  included:  tem¬ 
perature  and  salinity  (Sea-Bird  SBE21),  and  the  partial 
pressure  of  carbon  dioxide  (pC02;  General  Oceanics).  The 
ship’s  CTD  system  (Sea-Bird  9/11+)  included  sensors  for 
measuring  pressure,  temperature  (dual  sensors),  conductivity 
(dual  sensors),  and  dissolved  oxygen  (dual  sensors)  at  24 
Hz.  Data  were  processed  into  1  m  vertical  bins.  The  CTD 
data  were  used  on  site  for  identifying  features  including  the 
depths  of  the  mixed  layer,  deep  chlorophyll  maximum  layer, 
and  particle  maximum  layer. 

[21]  The  physical  context  of  the  specialized  optical  pack¬ 
age  called  the  MASCOT  (Multi-Angle  Scattering  Optical 
Tool),  which  is  described  below,  was  quantified  using  an 
SBE49  CTD  that  was  integrated  into  the  MASCOT  package 
[Twardowski  et  al.,  2012].  The  physical  data  collected  from 
the  R/V  Kilo  Moana’s  CTD  and  the  MASCOT  were  con¬ 
sistent  with  each  other.  Thus,  some  of  the  R/V  Kilo  Moana 
and  MASCOT  profiles  have  been  synthesized  for  plots 
shown  in  this  paper.  In  addition,  a  CTD  (Applied  Micro¬ 
systems)  was  deployed  as  part  of  the  HyperPro  optical  pro¬ 
filing  system  by  Lewis  et  al.  [201 1], 

[22]  One  of  the  RaDyO  program’s  goals  was  to  resolve 
and  model  the  dynamics  of  bubbles,  their  size  distributions, 
and  their  effects  on  the  subsurface  optical  fields.  For  this 
application,  an  acoustic  resonator  was  attached  for  some  of 
the  MASCOT  deployments  to  obtain  simultaneous  optical 
and  acoustic  data  to  ascertain  the  optical  and  acoustical 
effects  of  different  bubble  size  distributions  [Farmer  et  al., 
1998;  Czerski  et  al.,  2011;  Vagle  et  al.,  2012;  Twardowski 
et  al.,  2012].  One  of  the  challenges  of  this  work  was  to 
extend  the  frequency  range  of  the  acoustic  resonator  to  be 
able  to  measure  bubbles  of  radius  less  than  about  20  microns 
and  to  invert  attenuation  data  to  obtain  bubble  numbers.  In 
addition,  the  effort  included  a  goal  of  simultaneously  mea¬ 
suring  bubbles  with  the  acoustic  resonator  and  an  optical 
volume  scattering  instrument  [Czerski  et  al.,  2011;  Vagle 


et  al.,  2012;  Twardowski  et  al.,  2012]  to  determine  their 
correspondence  in  overlapping  bubble  size  ranges. 

[23]  Vertical  profiles  of  currents  were  provided  by  the  R/V 
Kilo  Moana’s  38  kHz  and  300  kHz  ADCP  systems  (RDI) 
whose  transducers  were  mounted  at  8  m  depth.  The  38  kHz 
system  obtained  data  downward  from  about  10  m  for  our 
study  and  the  300  kHz  system  obtained  data  from  about 
1 0  m  to  about  1 00  m.  Data  were  recorded  as  5  min  averages 
and  vertical  binning  of  data  was  accomplished  after  the 
cruise.  A  scanning  laser  wave  height  (lidar)  measurement 
system  was  deployed  from  a  boom  connected  between  the 
bows  of  the  R/V  Kilo  Moana’s  twin  hulls  [Zappa  et  al., 
2012],  The  lidar  altimeter  was  used  to  measure  large-scale 
wave  geometry  (height  and  slope)  by  scanning  downward 
toward  the  sea  surface  at  75  Hz  from  ^9  m  above  the  mean 
water  level  resulting  in  a  sampling  footprint  of  about  1 0  cm 
diameter  with  height  accuracy  of  about  ±5  mm  (rms). 

2.3. 1.3.  Optical  Measurements 

[24]  The  R/V  Kilo  Moana’s  underway  sampling  system 
measured  chlorophyll  a  fluorescence  (Turner  10-AU-005)  as 
part  of  the  data  stream  for  the  physical  parameters.  The 
ship’s  CTD  system  was  also  used  to  measure  optical  beam 
transmission  (WET  Labs  C-Star),  optical  backscatter  (WET 
Labs  LSS  6000),  chlorophyll  a  fluorescence  (Seapoint),  and 
PAR  (Biospherical  Instruments)  at  24  Hz  with  binning  at 
1  m  vertical  intervals.  An  optical  scattering  instrument 
(LISST  100X,  Sequoia)  was  also  mounted  on  the  CTD  to 
measure  near  forward  scattering  with  32  ring  detectors  within 
the  angular  range  of  0.08  to  13.5°  at  a  wavelength  of  532  nm. 
This  instrument  was  also  used  in  bench  top  mode  to  analyze 
discrete  water  samples.  Water  samples  were  collected  using 
the  CTD’s  24-bottle  rosette  for  analyses  of  absorption  and 
scattering  properties,  particle  size  distributions,  and  the  con¬ 
centrations  of  particulate  organic  carbon,  chlorophyll  a,  and 
total  dry  mass  of  suspended  particulate  matter.  A  Point- 
Source  Integrating  Cavity  Absorption  Meter,  multiple  path 
length,  liquid  core  waveguide  (MPLCW)  absorption  system, 
LISST  100X,  and  Coulter  counter  (Multisizer  III,  Beckman- 
Coulter)  were  used  for  some  of  this  work.  Similar  analyses 
were  made  using  water  from  the  surface  microlayer  collected 
from  a  glass  plate  using  a  small  boat  launched  from  the  R/V 
Kilo  Moana  [Wurl  et  al.,  201 1], 

[25]  One  of  the  objectives  of  RaDyO  was  to  obtain  high- 
fidelity  inherent  optical  property  (IOP)  data,  including  the 
volume  scattering  function,  in  order  to  compute  time- 
dependent  radiance  fields  from  incident  light  fields,  which 
are  modulated  by  the  surface  topography  and  its  transmis¬ 
sive  characteristics,  using  IOP  data  and  radiative  transfer 
models.  To  satisfy  this  objective,  Twardowski  et  al.  [2012] 
developed  a  special  optical  sampling  system  called  MAS¬ 
COT.  Briefly,  MASCOT  uses  a  658  nm  laser  light  source 
and  measures  the  volume  scattering  function  (VSF  or  j3)  at 
17  different  angles  from  10  to  170°  with  respect  to  the 
incident  beam  of  light  at  10°  intervals  with  a  sampling  rate 
of  20  Hz.  A  linear  polarizer  was  also  used  to  obtain  polarized 
VSF  measurements  in  order  to  increase  the  information 
content  concerning  the  particles’  characteristics  and  size 
distributions.  Other  scattering  measurements  were  obtained 
using  an  ECO-VSF  (WET  Labs)  instrument  for  angles  of 
100,  125,  and  150°  at  650  nm,  an  ECO-BB3  (WET  Labs) 
instrument  for  an  angle  of  1 17°  at  wavelengths  of  462,  530, 
and  657  nm  [see  Sullivan  et  al.  [2005]  and  Twardowski  et  al. 
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[2007b]  for  methodological  details),  and  a  LISST-100X 
Type  B  (Sequoia)  for  near  forward  scattering  in  32  log¬ 
spaced  angular  increments  from  0.08  to  13°  at  a  wavelength 
of  650  nm.  Light  absorption  coefficient,  a,  and  total  beam 
attenuation  coefficient,  c,  (with  total  scattering  coefficient,  b, 
computed  as  the  difference,  b  =  c  -  a)  measurements  were 
made  using  WET  Labs  ac-9  and  ac-s  instruments 
[Twardowski  et  al.,  1999],  These  instruments  sampled  at 
6  Hz,  with  the  former  measuring  at  9  wavelengths  and  the 
latter  at  87  wavelengths  in  the  visible. 

[26]  In  addition,  an  acoustic  resonator,  which  was 
described  earlier  [Czerski  et  al.,  2011;  Vagle  et  al.,  2012], 
and  a  radiance  camera  (RadCam,  described  next  and  by 
Lewis  et  al.  [2011])  were  deployed  with  MASCOT  for  some 
measurements.  Good  agreement  between  the  optical  and 
acoustic  attenuation  measurements  of  bubble  events  was 
obtained  [ Twardowski  et  al.,  2012;  Czerski  et  al.,  201 1], 

[27]  Lewis  et  al.  [2011]  deployed  new  instruments  capable 
of  measuring  the  full  radiance  distribution  above  and  below 
the  water  surface.  The  first  of  these,  the  HyperPro  (or 
HyperProbe,  Satlantic  Inc.),  is  a  profiling  device  that  is 
deployed  as  a  free-fall  vehicle.  The  HyperPro,  which  sends 
data  to  the  surface  for  real-time  data  processing,  measured 
the  following:  downwelling  (and  uplooking)  hyperspectral 
(124  wavelengths  from  350  to  700  nm  with  3  nm  between 
bands)  planar  irradiance,  upwelling  (and  downlooking) 
hyperspectral  radiance  with  spectral  characteristics  as 
described  for  the  downwelling  sensors  and  with  a  half¬ 
maximum  acceptance  angle  of  10°  (half-angle  of  Field- 
of-View).  A  CTD  (Applied  Microsystems),  chlorophyll  a 
fluorometer  and  CDOM  fluorometer  (WET  Labs),  2-channel 
backscatter  sensors  (470  and  670  nm),  a  GPS  position  and 
time  system,  and  2-axis  tilt  and  roll  sensors  were  also 
included  in  the  HyperPro  instrument  package.  This  system 
was  profiled  to  a  depth  of  60  m  and  also  operated  at  a  fixed 
depth  of  10  m  to  obtain  different  spatial  and  temporal  sta¬ 
tistics.  A  Tethered  Spectro-Radiometer  Buoy  (TSRB), 
which  is  essentially  the  same  measurement  device  as  the 
HyperPro,  was  tethered  to  a  surface  buoy  [Lewis  et  al, 
2011].  The  difference  being  that  the  downwelling  irradi¬ 
ance  sensor  is  inverted  to  measure  upwelling  irradiance.  This 
system  was  floated  away  from  the  ship  as  data  were  col¬ 
lected  to  avoid  ship  shadowing  and  reflecting  effects. 

[28]  The  Radiance  Camera  (RadCam,  Satlantic  Inc.)  is  a 
new  camera  designed  to  measure  spatial  radiance  distribu¬ 
tions  both  at  the  ocean  surface  from  the  ship’s  deck  and  at 
depth  [Lewis  et  al.,  2011],  A  surface  RadCam  instrument 
was  used  as  a  reference  camera  for  measuring  incident 
downwelling  radiance  distributions  (sky  radiance).  The 
subsurface  RadCam,  whose  characteristics  match  those  of 
the  deck-mounted  instrument,  includes  a  fish-eye  lens  and 
measured  downwelling  radiance  distributions  as  a  function 
of  depth,  time,  and  wavelength  or  Ld(z,  9,  <j).  A,  t)  and  Lu(z, 
9,  <f>.  A,  t)  where  z  is  depth,  9  is  the  zenith  angle,  <j>  is  the 
azimuthal  angle,  A  is  the  wavelength  of  light,  and  t  is  time. 
The  RadCam  measured  radiation  at  555  nm  (20  nm  band¬ 
pass  filter)  over  upward  looking  and  downward  looking 
hemispheres  with  an  angular  resolution  of  less  than  1°, 
sampling  rate  of  7  Hz,  and  a  dynamic  measurement  range  of 
106.  These  data  are  used  to  compute  (via  integration  over 
appropriate  solid  angles)  planar  irradiance,  scalar  irradiance, 
reflectance,  and  average  cosine,  and  inherent  optical 


properties  including  the  absorption  coefficient  and  scattering 
phase  function. 

[29]  The  RadCam  used  for  in- water  profiling  was  a  teth¬ 
ered  free-fall  system,  designed  to  minimize  the  perturbing 
effects  of  ship  shadowing  and  reflections.  The  profiling 
system  also  included  sensors  for  measuring  downwelling 
planar  irradiance,  Ed,  and  upwelling  nadir  radiance,  L u(9  = 
180°),  both  at  wavelengths  of  412,  443,  490,  and  555  nm, 
beam  attenuation  coefficient  (25  cm  path  length  at  532  nm 
(WET  Labs),  a  CTD  (Falmouth  Scientific),  and  a  tilt-heading 
sensor.  This  system  was  profiled  and  operated  at  fixed  depths 
of  2  to7  m.  Another  self-contained  RadCam  was  deployed 
from  the  Odyssey  AUV  described  later. 

[30]  The  R/V  Kilo  Moana  was  also  used  to  launch  a  sur¬ 
face  skimming  sampler  {Lil  KM)  and  two  autonomous 
underwater  vehicles  (AUVs).  These  platforms  are  described 
along  with  their  measurement  systems  in  the  following 
subsections. 

2.3.2.  Lil  KM 

[31]  Small-scale,  near-surface  physical,  chemical,  and 
optical  data  were  collected  using  a  specially  designed 
microlayer  skimmer,  dubbed  Lil  KM  (Figure  1),  during  the 
SBC  experiment  only  and  with  the  R/V  Kilo  Moana’’ s  small 
boat  during  both  experiments  [Wurl  et  al.,  2011],  The 
microlayer  skimmer  is  a  self-contained  remotely  operated 
double-hulled  vessel  1.5  m  long  and  1.4  m  wide  and  with  a 
draft  of  approximately  0.1  m.  The  skimmer  uses  10  rotating 
glass  disks  equipped  with  Teflon  wipers,  a  collection  mani¬ 
fold,  and  a  peristaltic  pump  to  collect  microlayer  samples 
from  eight  pairs  of  individually  controllable  250  ml  plastic 
bottles  for  further  analyses  in  the  laboratory.  The  duration  of 
each  mission  was  typically  between  1  and  2  h. 

2.3.2. 1.  Physical  Measurements 

[32]  During  the  SBC  experiment,  the  Lil  KM  skimmer  was 
equipped  with  an  additional  1.7  m  vertical  ‘mast’  suspended 
below  the  two  hulls  to  make  measurements  of  temperature  at 
2-5  depths  between  0.05  m  and  1.68  m  below  the  mean 
surface. 

2.3.2.2.  Optical  Measurements 

[33]  The  skimmer  also  measured  chlorophyll  a,  colored 
dissolved  organic  material  or  chromophoric  dissolved 
organic  matter  (CDOM),  and  660  nm  scattering  at  1.68  m 
depth  using  WET  Labs  ECO  triplet  sensors  mounted  on  the 
subsurface  mast. 

2.3.2.3.  Chemical  Measurements 

[34]  Surface-active  substances  modulate  the  microlayer  of 
the  near-surface  light  field.  Importantly,  compounds  con¬ 
tributing  to  CDOM  are  often  enriched  in  the  sea-surface 
microlayer  (SML).  These  substances  absorb  ultraviolet  and 
other  short-wavelength  light  energy  [ Wurl  et  al.,  2009]  and 
are  thus  important  for  RaDyO.  Lil  KM,  which  was  used 
for  microlayer  sampling  as  well  as  other  measurements, 
was  deployed  during  the  SBC  experiment  six  times  from 
September  14-22,  2008  for  periods  of  about  60-90  min, 
typically  around  1530  UTC  [Wurl  et  al.,  2011],  Microlayer 
and  bulk  water  samples  were  also  collected  from  the  R/V 
Kilo  Moana’s  small  boat  using  a  glass  plate  and  a  small 
peristaltic  pump  [Wurl  et  al.,  2009],  and  from  the  R/V  Kilo 
Moana’s,  CTD/rosette  system.  Some  of  the  SBC  data  were 
collected  within  visible  surface  slicks  and  near  R/P  FLIP 
(i.e.,  September  21,  2008)  over  the  2  week  sampling  period. 
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The  samples  were  analyzed  either  onboard  the  ship  or  taken 
back  to  a  land-based  laboratory  for  analyses. 

[35]  Particle  size  distributions  and  optical  absoiption  and 
scattering  by  particles  were  measured  to  characterize  particle 
properties  within  the  microlayer  in  collaboration  with 
Stramski  and  coworkers.  Analyses  were  also  done  for  dis¬ 
solved  carbohydrates,  surface-active  substances  (SAS,  e.g., 
surfactants),  and  surface  active  gel-particles  (transparent 
exopolymers,  TEP)  [Wurl  et  al.,  2009].  Wurl  et  al.  [2011] 
also  studied  the  removal  of  CDOM  from  the  SML  via  pho¬ 
tochemical  transformation  and  degradation  as  well  as  its 
enrichment  in  the  SML  via  rising  bubbles  which  scavenge 
organic  material  and  carry  it  into  the  SML.  Because  of  more 
adverse  wind  and  wave  conditions  off  Hawaii,  Lil  KM  was 
not  deployed.  However,  boat-based  measurements  similar  to 
those  done  during  the  SBC  experiment  were  made. 

2.3.3.  Odyssey  Autonomous  Underwater  Vehicle 

[36]  An  Odyssey  autonomous  underwater  vehicle  (AUV; 
hereafter  called  Odyssey )  [Moline  et  al,  2012,  Figure  1]  was 
deployed  with  physical  and  optical  instrumentation  during 
the  SBC  experiment  to  obtain  spatial  variability  data  in  the 
vicinity  of  the  R/P  FLIP  and  the  R/V  Kilo  Moana.  One  of 
the  advantages  of  AUV’s  is  that  they  can  collect  data  away 
from  ship’s  influence  (e.g.,  shading  and  reflections). 
Because  of  anticipated  adverse  wind  and  sea  states,  the 
Odyssey  AUV  was  not  deployed  during  the  Hawaii  experi¬ 
ment.  The  vehicle  is  a  free-flood,  modular,  midsized  plat¬ 
form  that  is  0.5  m  in  diameter  and  3.2  m  long  as  configured 
for  this  experiment.  The  Odyssey  can  sample  for  10  h  with 
operating  speeds  between  ~1.5  and  2  m/sec.  Pitch  and  roll 
near  the  surface  is  about  2°. 

[37]  During  the  SBC  experiment,  the  Odyssey  AUV  was 
operated  in  saw-tooth  patterns  with  10°  pitch  and  level 
(along  isobars)  flights  as  well  as  vertical  profiling  mode.  A 
10°  pitch  can  produce  quasi-vertical  profiles  in  the  upper 
30  m  for  every  170  m  horizontal  distance.  Level  flights 
provide  horizontal  spatial  scales  of  physical  and  optical 
fields,  while  saw-tooth  or  profiling  flights  also  provide  ver¬ 
tical  scale  and  gradient  information. 

2.3.3. 1.  Physical  Measurements 

[38]  The  Odyssey  was  equipped  with  a  pumped  CTD 
instrument,  downward  looking  300  kHz  ADCP,  micro¬ 
structure  sensing  package,  and  optical  sensors.  Turbulence 
measurements  were  made  with  a  micro  structure  package, 
which  consisted  of  two  fast  response  temperature  sensors, 
two  shear  probes,  pressure  sensor,  and  three-axis  linear 
accelerometer.  High  sampling  rates  were  needed  to  achieve 
small  spatial-scale  irradiance  measurements.  In  particular, 
turbulence  variables  and  irradiance  were  sampled  at  512  Hz 
and  pressure  was  sample  at  32  Hz. 

2.3.3.2.  Optical  Measurements 

[39]  A  nine-channel  absorption  and  beam  attenuation 
meter  (WET  Labs  AC-9)  was  mounted  in  the  midbody 
payload  section  of  the  Odyssey,  which  sampled  water 
pumped  from  the  same  location  as  the  CTD.  An  upward 
looking  7-channel  irradiance  sensor,  two  irradiance  sensors, 
and  a  downward  looking  optical  scattering  sensor  were  also 
mounted  in  the  payload  section.  To  achieve  high-resolution 
irradiance  measurements,  two  single  wavelength  irradiance 
sensors  (490  nm  and  532  nm)  were  integrated  into  the  RSI 
microstructure  package.  Sampling  rates  of  the  AC-9,  7- 
channel  irradiance  sensor,  and  backscattering  sensor  were 


8  Hz.  In  addition,  a  RadCam  (described  above)  was  mounted 
on  the  top  of  the  Odyssey  to  collect  radiance  data.  By  flying 
the  AUV  along  isobars,  it  was  possible  to  compute  power 
spectra  and  other  statistics  of  irradiance  fluctuations. 

2.3.4.  REMUS  Autonomous  Underwater  Vehicle 

[40]  A  REMUS  AUV  [Moline  et  al.,  2012]  (Figure  1)  was 
deployed  with  physical  and  optical  instrumentation  during 
both  the  SBC  and  Hawaii  experiments  to  obtain  spatial  data 
in  the  vicinity  of  the  R/P  FLIP  and  the  R/V  Kilo  Moana.  The 
vehicle  uses  a  single  pressure  hull  modular  design  and  is 
0.2  m  in  diameter  and  approximately  2  m  long.  The  vehicle 
can  sample  at  1.5  m/sec  for  10  h  in  the  temperature  envi¬ 
ronment  of  the  study  sites. 

[41]  During  the  SBC  and  Hawaii  experiments,  the  REMUS 
AUV  was  operated  in  both  fixed-level  stepwise  flight  and 
saw-tooth  mode  with  a  similar  pitch  as  the  Odyssey.  Level 
flights  provided  horizontal  spatial  scales  of  physical  and 
optical  fields,  while  saw-tooth  or  profiling  flights  also  pro¬ 
vided  vertical  scales  and  gradients. 

2.3.4. 1.  Physics 

[42]  The  REMUS  AUV  was  equipped  with  a  CTD  and 
upward  and  downward  looking  1200  kHz  ADCPs.  The  CTD 
sensors,  which  were  mounted  in  the  nose  of  the  vehicle, 
sampled  at  12  Hz  and  were  able  to  sample  vertical  gradients 
in  the  physical  parameters  at  resolutions  better  than  2  cm. 

2.3.4.2.  Optics 

[43]  Also  mounted  in  the  nose  section  of  the  REMUS 
AUV  with  the  CTD  sensor  were  a  7-channel  downward 
irradiance  meter  and  sensors  for  measuring  upward  radiance, 
chlorophyll  a,  optical  backscatter,  and  CDOM. 

2.3.5.  R/P  FLIP 

[44]  For  the  SBC  experiment,  R/P  FLIP  was  used  for  data 
collection  from  September  1 1  until  September  22,  2008  at 
34.2053°  N,  119.6288°  °W.  For  the  Hawaii  experiment 
south  and  southwest  of  the  island  of  Hawaii,  the  free  drifting 
R/P  FLIP  sampled  from  September  1  until  September  14, 
2009. 

2.3.5. 1.  Meteorological  Measurements 

[45]  Two  separate,  but  similar,  sets  of  meteorological  data 
were  collected  during  both  of  the  RaDyO  field  experiments 
by  Zappa  et  al.  [2012]  and  Melville  et  al.  (manuscript  in 
preparation,  2012).  Measurements  included:  near-surface 
barometric  pressure,  wind  speed  and  direction,  relative  and 
specific  humidity,  air  temperature,  fast  response  water 
vapor,  C02,  and  temperature  sensors  for  heat  flux  measure¬ 
ments,  longwave  radiative  downwelling  flux,  and  shortwave 
downwelling  irradiance  (Kipp  and  Zonen  pyrgeometer).  All 
sensors  were  mounted  at  10  m  above  sea  level.  Postproces¬ 
sing  of  data  allowed  determinations  of  momentum  flux 
(wind  stress),  turbulent  latent  and  sensible  heat  fluxes,  and 
net  heat  flux. 

2.3. 5.2.  Physical  Measurements 

[46]  R/P  FLIP  was  used  by  the  UCSB  group  for  several 
physical  and  optical  profiling  and  fixed-depth  measure¬ 
ments.  During  the  SBC  experiment,  a  UCSB  optical-physi¬ 
cal  profiling  system  was  deployed  to  depths  of  15  or  30  m. 
The  profiling  system  included  a  CTD  along  with  several 
optical  sensors  described  in  the  next  section.  The  physical 
and  optical  data  were  depth-binned  to  10  cm  and  time- 
binned  to  1  Hz.  A  total  of  30  casts  were  made  during  the 
field  experiment:  15  profiles  down  to  a  depth  of  30  m,  and 
15  10  min  long  time  series  at  a  depth  of  2  m.  During  the 
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vertical  profiling,  the  package  was  left  at  30  m  for  5  min  and 
was  brought  back  to  the  surface  at  a  rate  of  1  m/min.  Mea¬ 
surements  were  done  from  September  1 1  until  September 
14.  Two  vertical  profiles  were  taken  each  day,  one  between 
1200  and  1300  and  another  between  1530  and  1630; 
whereas  the  time  series  at  2  m  were  collected  every  hour, 
usually  from  900  to  1700.  Deep  profiles  were  made  at  1200, 
1330,  and  1500  local  time  on  September  17  and  at  1000, 
1200,  1400,  1600,  and  1800  on  September  20. 

[47]  No  profile  measurements  were  made  with  the  UCSB 
profiling  system  during  the  Hawaii  experiment.  However, 
several  UCSB  physical  and  optical  sensors  were  mounted  on 
R/P  FLIP’S  hull.  Specifically,  22  temperature,  4  salinity,  and 
2  pressure  sensors  made  measurements  in  the  upper  85  m 
with  these  instruments  (Figure  2). 

[48]  Vertical  profiles  of  currents  and  wave  directional 
spectra  were  obtained  by  a  600  kHz  wave-enabled  upward 
looking  600  kHz  RDI ADCP  mounted  on  R/P  FLIP’S  hull  at 
30m  water  depth  during  the  SBC  experiment,  and  wave- 
enabled  upward  looking  600  kHz  and  downward  looking 
300  kHz  RDI  ADCPs  mounted  at  14  m  water  depth  during 
the  Hawaii  experiment  (Figure  2)  were  deployed  by  Melville 
et  al.  (manuscript  in  preparation,  2012). 

[49]  Vagle  et  al  [2012]  deployed  a  CTD  package  mounted 
at  30  m  on  R/P  FLIP  to  collect  data  every  10  min  throughout 
the  experimental  periods.  The  package  also  included  a  dis¬ 
solved  oxygen  sensor  and  a  gas  tension  device.  These  data 
were  used  for  interpretation  of  their  bubble  measurements  as 
well  as  upper  ocean  physics.  In  addition,  hull-mounted 
internally  recording  thermistors  were  deployed  at  depths  of 
3,  7,  15,  31,  and  63  m. 

[50]  Surface  roughness  variability,  small-scale  waves,  and 
wave  breaking  events  were  observed  using  suites  of  instru¬ 
ments  measuring  on  scales  down  to  millimeters  from  R/P 
FLIP  during  both  field  experiments  by  two  wave-surface 
roughness  groups:  Zappa  et  al.  [2012]  and  Melville  et  al. 
(manuscript  in  preparation,  2012).  These  data  sets  were 
collected  to  examine:  (1)  wave  height  and  slope,  (2)  small- 
scale  sea-slope  field  topography,  and  (3)  microbreaker  and 
whitecap  crest  length  spectral  density  of  breaker  propagation 
speed  as  functions  of  wind  speed,  wind  stress,  and  dominant 
sea  states. 

[51]  Zappa  et  al.  [2012]  deployed  a  moderate  field-of- 
view  CED1P  infrared  (IR)  camera,  three  digital  video  cam¬ 
eras,  and  a  Riegl  laser  altimeter  as  a  system  from  R/P  FLIP’ s 
starboard  boom.  The  IR  camera  measured  thermal  radiation 
at  wavelengths  of  7. 7-9. 3  /mi  emitted  by  the  ocean  surface 
using  a  320  by  240  MCT  focal -plane  array  with  14-bit  dig¬ 
itization  sampled  at  60  Hz.  Resulting  temperature  resolution 
was  0.02°C  and  calibration  was  better  than  0.05°C.  The 
three  digital  visible  CCD  cameras  recorded  images  sampling 
at  20  Hz,  two  at  1250  x  1360  resolution  and  12-bit  digiti¬ 
zation  depth,  and  one  at  1000  x  1000  resolution  with  12-bit 
digitization  depth.  The  laser  altimeter,  which  operated  at  a 
wavelength  of  0.905  pm,  measured  the  distance  to  the  sur¬ 
face.  The  IR  and  visible  camera  images  were  postprocessed 
to  compute  statistics  on  the  scale,  frequency,  and  speed  of 
microbreaking  and  whitecapping  events  from  scales  of  order 
0.1  m  s_1  up  to  scales  of  order  10  ms-1.  In  particular,  these 
data  were  used  for  determining  breaking  crest  length  spectral 
density  distributions  and  their  higher  moments. 


[52]  The  laser  altimeter  data  were  used  to  obtain  wave 
spectra  and  other  wave  statistics  including  significant  wave 
height  and  wave  frequency.  In  addition,  a  synchronous 
orthogonal  75  Hz  linear  scanning  laser  altimeter  system  was 
used  to  provide  spatiotemporal  properties  (interlaced  x-z  and 
y-z  scans)  of  the  wave  height  field  resolved  to  the  order  of 
0.5  m  wavelengths.  These  measurements  also  provide  phase 
of  polarimeter  images  (described  below)  of  surface  micro¬ 
structure  needed  for  quantifying  short  wave  modulation. 
This  system  was  also  located  at  10  m  height.  The  system 
scanned  in  two  dimensions  (using  Cartesian  coordinates) 
and  the  sampling  footprint  was  about  10  cm.  The  accuracy 
of  the  measured  height  was  about  ±5  mm  (mis).  The  three 
cameras  were  used  to  observe  wave  breaking,  foam,  and 
whitecapping  over  fields  of  10  x  15  m  (2  cameras)  and 
100  x  200  m  ( 1  camera),  respectively.  The  former  cameras 
were  mounted  near  the  fixed  and  scanning  lidar  system  close 
to  the  end  of  a  boom  (10m  height)  to  measure  intermediate 
scale  breakers  while  the  other  was  mounted  on  R/P  FLIP’ s 
crow’s  nest  at  about  26  m  above  water  level  for  a  broader 
viewing  angle  to  record  larger-scale  breaking  events.  The 
sampling  rates  for  all  cameras  were  either  10  or  20  frames 
per  sec.  Data  were  generally  collected  for  20  min  every  hour; 
however,  data  were  recorded  more  often  during  periods  of 
more  frequent  wave  breaking  events. 

[53]  The  phase-resolved,  spatial-temporal  history  of  small 
waves  was  measured  using  the  shape-from-polarimetry 
(SFP)  technique  first  described  in  Zappa  et  al.  [2008].  The 
SFP  technique  relates  the  change  in  polarization  of  skylight 
reflecting  from  water  to  infer  the  orientation  of  the  water 
surface  at  the  point  of  reflection. 

[54]  With  funding  from  an  Office  of  Naval  Research 
DURIP  award,  and  in  cooperation  with  Polaris  Sensor 
Technologies,  an  imaging  polarimeter  was  specifically 
designed  and  built  for  oceanographic  applications  by  Zappa 
et  al.  [2012],  The  imaging  polarimeter  collected  782  x  582 
pixel  10-bit  monochrome  images  at  60  Hz.  The  instrument 
was  equipped  with  a  ~3.5  deg  field-of-view  lens,  which 
resulted  in  a  surface  footprint  of  approximately  lxl  m.  The 
design  and  results  are  described  in  more  detail  in  Zappa 
et  al.  [2012], 

[55]  The  imaging  polarimeter  was  deployed  in  a  water 
resistant  housing  from  the  starboard  boom  9  m  above  the 
surface.  The  sensor  was  oriented  so  that  it  observed  the 
ocean  surface  at  a  37  deg  incidence  angle.  The  housing 
contained  the  imaging  polarimeter,  an  attitude  heading  and 
reference  system  with  a  dual  GPS  antenna  for  collecting 
heading  data,  a  Camera  Link  to  fiber  optic  interface,  and  a 
temperature  sensor.  The  raw  noncompressed  image  data  and 
orientation  data  were  tagged  with  GPS  time  and  recorded  in 
real-time  using  an  IO  Industries  frame  grabber  and  SAS 
RAIS  system. 

[56]  Melville  et  al.  (manuscript  in  preparation,  2012) 
deployed  a  set  of  electro-optical  sensors  from  the  starboard 
(for  SBC)  and  port  (for  Hawaii)  booms  of  R/P  FLIP  to 
characterize  the  influence  of  wave  breaking  at  the  ocean 
surface  on  oceanic  radiance  and  imaging.  These  included  a 
Riegl  Q240i  scanning  laser  altimeter  for  surface  displace¬ 
ment  and  slope  measurements,  stereo  visible  and  infrared 
cameras  (JaiPulnix  TM4100CL  4Mpx  visible  cameras  and 
FL1R  SC6000  LWIR  infrared  cameras  (640  x  512px)  sam¬ 
pled  at  10  Hz  and  40  Hz  respectively)  for  3D  surface 
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reconstruction,  wave  kinematics  and  breaking  characteristics 
estimation,  dual  11  megapixel,  12  bit  digital  cameras 
mounted  on  the  crow’s  nest  of  R/P  FLIP  (Figure  2)  for 
whitecap  statistics  and  kinematics.  In  addition,  an  electro- 
optical  package  was  deployed  at  the  very  end  of  the  port 
boom  during  the  Hawaii  deployment.  [This  instrument  was 
not  deployed  during  the  SBC  experiment  after  the  failure  of 
the  port  boom,  subsequently  replaced  prior  to  the  Hawaii 
deployment.]  For  this  experiment  we  upgraded  the  midwave 
infrared  (MWIR)  digital  camera  to  an  Indigo  Phoenix  cam¬ 
era  (320  x  256px,  50  Hz  frame  rate),  and  replaced  the  video 
camera  by  a  JaiPulnix  3CCD  color  camera  (CV-M9CL, 
1360  x  1024px,  up  to  50  Hz  sampling  rate).  The  system  also 
includes  a  60-W  air-cooled  C02  laser  (Synrad  Firestar  T60) 
equipped  with  an  industrial  marking  head  (Synrad  FH  index) 
with  two  computer  controlled  scanning  mirrors  and  a  laser 
altimeter  (Riegl  LD90-3100-EHS).  To  characterize  the  light 
transmission  across  the  water  surface,  this  group  also 
deployed  a  downward  looking  LED  screen  (Figure  2)  used 
to  display  preprogrammed  targets,  imaged  through  the  water 
surface  by  an  upward  looking  color  camera  (Prosilica 
GC1380C)  mounted  on  a  small  underwater  profiler  fitted 
with  additional  subsurface  instruments  to  characterize  tur¬ 
bulence  and  bubble  injection  below  the  surface  (Nortek  HR 
profiler  and  Vector  ADV).  Instruments  were  sampling  con¬ 
tinuously,  except  the  scanning  lidar  and  imaging  systems 
that  were  recording  for  20  min  per  hour. 

[57]  Information  concerning  upper  ocean  bubble  structure 
was  obtained  by  Czerski  et  al.  [2011]  and  Vagle  et  al. 
[2012],  Data  were  obtained  within  a  radial  distance  of 
about  250  m  of  R/P  FLIP  using  a  Doppler  sonar  system 
consisting  of  four  100  kHz  side  scan  transducers  oriented 
orthogonally,  15°  up  from  the  horizontal  plane,  at  a  depth  of 
30  m.  A  horizontally  mounted  300  kHz  acoustic  Doppler 
current  profiler  (ADCP)  was  also  mounted  at  30  m  depth  and 
operated  throughout  the  experiment.  A  CTD  package 
described  earlier  was  mounted  at  30  m  to  collect  tempera¬ 
ture,  conductivity,  and  pressure  data  at  30  m  every  10  min.  A 
dissolved  oxygen  sensor  and  a  gas  tension  device  (GTD) 
were  included  in  the  CTD  package  for  measurements  of  the 
total  pressure  of  all  dissolved  gases.  Since  air  is  mostly  made 
up  of  the  biologically  active  oxygen  (02)  and  the  practically 
inert  nitrogen  (N2),  the  GTD  measurements  combined  with 
02  measurements  allow  for  estimating  N2  concentrations, 
which  are  needed  for  proper  interpretation  of  the  bubble 
measurements. 

[58]  For  in  situ  near-surface  measurements  of  optically 
active  properties  of  the  upper  ocean,  a  small  (1.5  m  by 
1.5  m)  wave-following  float  (Figure  2c)  was  also  deployed 
between  the  booms  of  R/P  FLIP  by  Vagle  et  al.  [2012].  The 
float  was  equipped  with  two  acoustical  resonators  for  in  situ 
measurements  of  the  bubble  size  distributions  at  two  depths 
(0.5  and  1.5  m),  three  2  MHz  Doppler  sonars  (Dopbeams) 
for  measuring  backscatter  and  near-surface  turbulence  along 
1  m  paths,  conductivity-temperature  packages  at  0.35  and 
1.8  m,  an  inductive  conductivity  sensor  used  to  measure 
high  air  fractions  in  breaking  waves,  and  a  sensor  package 
with  a  pressure  sensor,  plus  heading  and  tilt  sensors.  In 
addition,  Zappa  et  al.  [2012]  deployed  a  side-looking  high- 
resolution  Aquadopp,  which  was  mounted  at  2  m  depth  to 
measure  backscatter  and  near-surface  turbulence  along  1  m 
paths.  Attached  at  5  m  below  the  float  was  an  upward 


looking  200  kHz  short  pulse  (250  /js)  backscatter  sonar 
system  used  to  map  the  larger-scale  bubble  field  in  the 
vicinity  of  the  float.  Finally,  two  black  and  white  video 
cameras  collected  images  in  the  vicinity  of  the  float  to 
ascertain  the  presence  or  absence  of  wave  breaking  over  the 
float  and  to  indicate  the  presence  or  absence  of  surface  foam 
lines  associated  with  Langmuir  circulation. 

2.3. 5.3.  Optical  Measurements 

[59]  The  UCSB  profiling  CTD  system,  which  was  used 
only  for  the  SBC  experiment,  also  included  a  hyperspectral 
absorption-attenuation  instrument  for  measuring  optical 
absorption  and  attenuation,  a  and  c  (ac-s,  87  wavelengths), 
another  absorption-attenuation  instrument  (WET  Labs 
AC-9)  for  measuring  a  and  c  with  a  filter  on  the  a-tube  for 
CDOM  (Nucleopore  0.2  //,m  pore  size  filter,  9  wavelengths  in 
the  visible),  an  optical  backscatter  instrument  for  measuring 
bb  at  9  wavelengths  in  the  visible,  a  fluorometer-turbidity 
instrument  for  measuring  chlorophyll  fluorescence  (470  nm 
ex/695  nm  em)  and  turbidity  (700  nm),  and  a  near  forward 
scattering  instrument  for  measuring  particle  size  distribu¬ 
tions  (characteristics  similar  to  those  of  the  LIS  ST  instru¬ 
ments  described  earlier). 

[60]  Fixed-depth  time  series  measurements  were  also 
made  with  UCSB  systems  during  the  SBC  experiment  using 
another  measurement  system  that  was  mounted  on  the  hull 
of  R/P  FLIP  at  a  depth  of  2  m.  This  system,  called  the  WQM 
(WET  Labs  [Janzen  et  al,  2008]),  included  a  combined 
CTD/dissolved  oxygen  system,  and  a  combined  chlorophyll 
a  fluorometer  and  turbidity  sensor.  For  the  Hawaii  experi¬ 
ment,  UCSB  R/P  FLIP  hull-mounted  optical  measurements 
were  made  using  WQM  instruments  as  described  above  for 
fluorescence  and  turbidity  at  four  depths  (13,  20,  40,  and 
68.5  m). 

[61]  Another  suite  of  specialized  optical  measurements 
was  utilized  during  both  the  SBC  and  Hawaii  experiments. 
Radiometric  measurements  of  wave-induced  light  fluctua¬ 
tions  were  made  with  the  Underwater  Porcupine  Radiometer 
System  [Darecki  et  al.,  2011],  The  Porcupine  is  a  new 
instrument  developed  specially  for  the  RaDyO  project.  It  has 
a  unique  design  and  capability  to  sample  rapid  fluctuations 
in  both  radiance  and  irradiance  fields  at  1  kHz.  The  Porcu¬ 
pine  was  equipped  with  7  radiometric  sensors  to  measure 
downwelling  plane  irradiance,  Ed,  at  wavelengths  of  365, 
410,  443,  488,  532,  610,  and  670  nm.  The  cosine  collectors 
were  2.5  mm  in  diameter.  In  addition,  the  Porcupine 
included  16  radiance  sensors  to  measure  downwelling  radi¬ 
ance,  Ld,  at  532  nm  at  different  zenith  angles  within  two 
orthogonal  azimuthal  planes. 

[62]  Other  radiometric  measurements  were  made  with 
several  Ramses  (TriOS  GmbH)  hyperspectral  radiometers 
[Darecki  et  al.,  2011],  which  provided  time-averaged  light 
field  characteristics  (averaging  times  typically  from  about 
0.5  s  to  a  few  seconds).  The  measurements  with  these 
radiometers  provide,  for  example,  hyperspectral  average 
cosine  data  for  the  underwater  light  field,  an  apparent  optical 
property  characterizing  the  angular  distribution  of  the  light 
field.  The  Ramses  radiometers  attached  to  the  Porcupine 
measured  in- water  hyperspectral  (190  channels  within  the 
wavelength  range  320-950  nm;  3.3  nm  resolution)  down¬ 
welling  plane  irradiance,  Ed,  downwelling  scalar  irradiance, 
Eod,  downwelling  radiance  (0  deg,  zenith  radiance),  Ld,  and 
chlorophyll  fluorescence.  A  radiometer  mounted  on  a 
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surface  float  was  used  to  measure  hyperspectral  (320-950 
nm)  upwelling  radiance,  Lu,  just  below  the  sea  surface  (180 
deg;  nadir  radiance),  and  two  deck  radiometers  measured 
hyperspectral  Ed  and  Eod  incident  on  the  sea  surface. 

[63]  Bhandari  et  al.  [2011]  used  cameras  designed  for 
measurements  of  polarized  radiance  distributions  for 
deployments  from  R/P  FLIP  above  the  surface  (SKY  CAM) 
and  in  the  water  (DPOL,  Downwelling  POL  camera  system) 
for  both  the  SBC  and  Hawaii  RaDyO  experiments.  The 
DPOL  system  allows  measurements  of  one  hemisphere  of 
the  polarized  radiance  distribution  at  5  wavelengths  (442, 
488,  520,  550,  589  nm)  and  can  be  deployed  for  either 
downwelling  or  upwelling  radiance  distributions.  For  the 
downwelling  radiance  distribution,  the  DPOL  was  deployed 
below  several  of  the  wave  measurement  systems  whereas  for 
the  upwelling  radiance  distribution,  the  system  was 
deployed  using  floats  away  from  the  R/P  FLIP.  By  using  a 
combination  of  four  simultaneously  acquired  images,  the 
four  Stokes  parameters  I,  Q,  U,  and  V  can  be  acquired  at  a 
given  wavelength.  When  measuring  the  downwelling  radi¬ 
ance  distribution,  depths  were  varied  from  0.5  to  40  m  in  the 
SBC  and  1  to  70  m  during  the  Hawaii  experiment,  with 
sampling  at  each  depth  for  approximately  30  min. 

[64]  The  SKY  CAM  measured  the  surface  downwelling 
spectral  radiance  distribution  at  8  wavelengths  (442,  488, 
520,  550,  589,  650,  766,  and  855  nm).  This  system  is  a 
newer  version  of  an  older  SKY  CAM  system  [  Voss  and  Liu, 
1997].  This  system  was  deployed  from  an  R/P  FLIP  boom 
and  the  camera  was  placed  in  a  gyroscopically  stabilized 
mount  with  azimuthal  control.  This  system  typically  col¬ 
lected  data  at  each  wavelength  every  1 5  min  throughout  the 
day.  With  this  system,  3  sequential  images  are  combined  to 
obtain  three  of  the  polarization  Stokes  vectors  (I,  Q,  and  U), 
since  in  the  atmosphere  the  Stokes  vector  V  is  very  small. 
Aerosol  optical  depth  (AOD)  was  also  measured  simulta¬ 
neously  with  DPOL  and  SKY  CAM  observations  with  a 
handheld  Sun  photometer  (Microtops).  During  the  Hawaii 
experiment,  a  RADCAM  [Lewis  et  al.,  2011]  was  also 
deployed  from  R/P  FLIP. 

2.3.6.  National  Data  Buoy  Center  -  Buoy  46053 

[65]  During  the  SBC  experiment,  meteorological  and 
oceanographic  data  were  routinely  collected  by  NOAA’s 
National  Data  Buoy  Center  (NDBC)  from  NDBC  buoy 
46053.  The  3  m  disc  buoy  was  located  in  the  SBC  about 
22.2  km  (12  nmi)  southwest  of  Santa  Barbara  at  34.25  N 
1 19.84  W  in  water  of  depth  450  m  or  about  22  km  generally 
west  of  the  RaDyO  SBC  study  site  (Figure  4).  Data  of 
interest  for  our  experiment  were  collected  1  September  to  1 
October  2008. 

2.3. 6.1.  Meteorological  Measurements 

[66]  Time  series  data  collected  from  the  NDBC  buoy 
included  wind  speed  and  direction  at  5  m  height  above  the 
sea  surface,  wind  gust,  air  temperature,  barometric  pressure, 
and  pressure  tendency  at  sea  level. 

2.3. 6.2.  Physical  Measurements 

[67]  NDBC  buoy  46053  oceanographic  time  series  data 
included:  ocean  temperature  at  a  depth  of  0.6  m,  mean  wave 
height,  significant  wave  height,  swell  height,  direction  and 
period,  mean  wave  direction,  average  wave  period,  and 
dominant  wave  period  with  the  maximum  wave  energy. 
More  information  on  these  measurements  can  be  found  at 
http://www.  ndbc.  noaa.gov/station  _page.php?station=46053. 


2.3.7.  Shore-Based  High-Frequency  Surface  Current 
Observations 

[68]  Near-surface  currents  in  the  Santa  Barbara  Channel 
have  been  measured  since  1997  using  a  shore-based  array  of 
high-frequency  (HF)  radars  along  the  mainland  coast  and  on 
Santa  Cruz  Island.  Several  SBC  studies  have  used  these 
radars  previously  [e.g.,  Nishimoto  and  Washburn,  2002; 
Beckenbach  and  Washburn,  2004;  Emeiy  et  al.,  2004; 
Bassin  et  al,  2005;  Anderson  et  al.,  2006,  2008],  The  HF 
radars  used  in  this  study  operated  at  12-13  MHz  and  mea¬ 
sured  currents  over  the  upper  1  m  of  the  water  column.  The 
radial  components  of  current  velocity  were  measured  every 
10  min  and  spatial  resolution  of  radial  currents  was  1.5  km 
in  range  and  5°  in  azimuth.  Averages  of  1  h  were  used  for 
the  plots  shown  in  this  paper.  Daily  plots  were  computed  as 
25  h  averages  to  minimize  effects  of  tidal  flows.  Surface 
current  vectors  were  interpolated  onto  a  2  km  square  grid 
based  on  all  radial  currents  within  3  km  of  each  grid  point 
using  a  least  square  fit.  An  eddy  tracking  algorithm  devel¬ 
oped  by  Nencioli  et  al.  [2010]  was  applied  to  the  HF  radar 
data  collected  during  the  experiment  and  results  are 
described  by  Dong  et  al.  [2012], 

2.3.8.  Aircraft  Operations 

[69]  The  airborne  scanning  L1DAR  system  developed  by 
Reineman  et  al.  [2009]  was  installed  in  the  cabin  of  a  Twin 
Piper  Comanche  aircraft  to  measure  the  spatial  evolution  of 
the  wavefield  in  the  vicinity  of  R/P  FLIP  and  R/V  Kilo 
Moana  during  the  Santa  Barbara  Channel  Experiment  on 
16  and  19  September  2008  (Melville  et  al.,  manuscript  in 
preparation,  2012). 

2.3.9.  Satellites 

[70]  In  order  to  characterize  the  general  setting  and  larger- 
scale  variability  of  the  SBC  and  Hawaii  RaDyO  experi¬ 
mental  regions,  satellite-based  wind,  sea  surface  tempera¬ 
ture,  and  ocean  color  were  obtained.  Surface  roughness 
(synthetic  aperture  radar  (SAR))  measurements  were  col¬ 
lected  for  the  SBC  experiment  only.  Satellite-derived  SST 
and  chlorophyll  a  imagery  were  at  times  not  useful  due  to 
cloud  cover. 

2.3.9. 1.  Wind  Measurements 

[71]  For  regional  surface  wind  speeds  and  directions, 
NASA’s  QuikSCAT  scatterometer  was  utilized.  QuikSCAT, 
which  is  a  polar  orbiting  satellite,  provided  wind  data  over 
an  1800  km  wide  swath  for  our  two  study  regions  (see 
QuikSCAT  website:  www.podaac.jpl.nasa.gov/quikscat).  The 
retrievals  of  wind  speed  and  direction  from  QuikSCAT  gave 
twice-daily  data  with  spatial  resolution  of  25  km  x  25  km  on 
the  Earth’s  surface. 

2.3.9.2.  Sea  Surface  Temperature  Measurements 

[72]  NASA’s  MOD1S  sea  surface  temperature  (SST) 
imagery  data  (9  km  pixel  resolution)  were  obtained.  The 
MODIS  Aqua  satellite  images  the  full  earth  every  1-2  days, 
therefore  the  daily  composites  usually  provided  coverage 
over  our  study  areas  every  other  day.  Eight  day  averages 
over  9x9  km  areas  were  used  for  our  analyses.  The  pre¬ 
ceding  data  were  level  3  products  obtained  from  podaac.jpl. 
nasa.gov. 

2.3.9.3.  Color  Measurements 

[73]  NASA’s  MODIS  ocean  color  data  were  used  to 
determine  average  near-surface  chlorophyll  a  concentrations 
over  8  days  and  areas  of  9  x  9  km  (data  obtained  from 
oceans.gsfc.nasa.gov).  The  MODIS  imagery  had  9  km 
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Figure  3.  Satellite  SAR  image  from  Envisat  ASAR  (copyright  ESA)  obtained  20  September  2008,  at 
05:45  UTC.  At  the  time  of  acquisition,  the  nearby  NOAA  buoy  46053  measured  a  wind  speed  of  9.4 
m/s  with  a  direction  of  280°,  wave  height  of  1.3  m,  wave  direction  of  284°  and  dominant  wave  period 
of  13.8  s.  The  high  winds  account  for  the  overall  bright/high  radar  backscatter  returns.  Swell  with  refrac¬ 
tion  patterns  around  Santa  Cruz  Island  are  detectable.  Atmospheric  frontal  patterns  are  discernible  at  the 
east  end  of  Santa  Cruz  and  in  the  northerly  midchannel.  Closer  to  the  coast,  natural  seeps  appear  as  dark, 
thin,  curvilinear  features,  including  a  pair  of  seeps  emanating  near  a  pair  of  drilling  platforms.  Dark  wind 
shadows  are  seen  close  to  the  coastline  from  west  of  Santa  Barbara  southward  toward  Ventura. 


ground  pixel  resolution  with  chlorophyll  a  concentrations 
being  calculated  using  standard  NASA  data  processing 
routines  and  algorithms  (see  http://oceancolor.gsfc.nasa.gov/ 
PRODUCTS/). 

2.3.9. 4.  Surface  Roughness-Synthetic  Aperture  Radar 
Measurements 

[74]  For  the  SBC  experiment,  satellite  synthetic  aperture 
radar  (SAR)  imagery  was  obtained  by  Ben  Holt  (JPL)  from 
the  European  Space  Agency’s  Envisat  and  ERS-2  satellites 
with  ASAR  and  SAR  sensors,  respectively.  Envisat  ASAR 
data  were  primarily  obtained  with  each  beam  having  25  m 
pixel  resolution  and  100  km  swath  width.  The  ERS-2  SAR 
provided  data  with  25  m  resolution  and  100  km  swath  width. 
Imagery  was  obtained  on  the  following  dates:  4,  7,  10,  20, 
23,  and  26  September  in  2008.  These  data  sets  are  useful  for 
detecting  surface  features  such  as  slicks  and  estimating  wave 
direction  and  wavelength,  wind  speed,  and  current  patterns 
[DiGiacomo  and  Holt,  2001:  DiGiacomo  et  al.,  2004].  A 
sample  SAR  image  frame  is  shown  in  Figure  3. 

3.  Background  for  the  Santa  Barbara  Channel 
RaDyO  Study 

[75]  The  Santa  Barbara  Channel  (SBC)  is  approximately 
40  by  100  km  in  size  and  reaches  a  maximum  depth  of  about 


600  m  near  the  center  of  the  basin  (bathymetry  shown  in 
Figure  4).  The  SBC  is  bordered  by  the  Santa  Ynez  Mountain 
Range  (reaching  over  1200  m  in  elevation)  to  the  north  and 
the  Santa  Barbara  Channel  Islands  to  the  south.  The  SBC  is  a 
relatively  well-defined  and  delimited  oceanic  water  body  as 
shown  in  Figure  4.  Local  topography  and  coastline  geometry 
significantly  influence  the  local  meteorological  conditions 
[e.g.,  Carvalho  et  al,  2012], 

[76]  The  Radiance  in  a  Dynamic  Ocean  (RaDyO)  program 
conducted  its  second  field  experiment  in  the  SBC  (Figure  4). 
Observational  platforms  including  the  R/V  Kilo  Moana,  R/P 
FLIP,  two  autonomous  underwater  vehicles  (AUVs),  a  sur¬ 
factant  skimmer,  and  a  small  aircraft  (see  Figure  1  for  all 
except  the  aircraft)  were  used  to  collect  data  during  the  period 
of  3-25  September  2008  in  the  vicinity  of  34.2053°  N, 
119.6288°  in  waters  of  depth  of  about  170  m.  These  plat¬ 
forms  collectively  obtained  high  temporal  and  spatial  reso¬ 
lution  data  sets.  During  the  RaDyO  Santa  Barbara  Channel 
experiment,  the  R/V  Kilo  Moana  operated  within  about 
2  km  (generally  north)  of  R/P  FLIP  (Figure  4)  during  day¬ 
light  hours,  but  had  to  leave  the  area  during  a  few  hours  in 
the  evenings.  R/P  FLIP  began  sampling  on  10  September 
2008  and  finished  on  25  September  2008  whereas  R/V  Kilo 
Moana  began  measurements  on  9  September  2008  and 
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Figure  4.  Geographic  and  bathymetric  map  of  region  of  the  RaDyO  field  experiment  conducted  in  Sep¬ 
tember  2008  in  the  Santa  Barbara  Channel.  The  main  site  of  the  RaDyO  SBC  experiment  is  indicated  by 
the  red  box  and  complementary  measurement  and  instrumentation  locations  are  described  in  the  key. 


terminated  them  on  22  September  2008.  Again,  comple¬ 
mentary  in  situ,  aircraft,  satellite  and  shore-based  data  sets 
have  been  analyzed  to  allow  broader-scale  interpretation  of 
the  dedicated  RaDyO  data  sets. 

3.1.  Results  for  the  Santa  Barbara  Channel  Experiment 

[77]  General  atmospheric  and  oceanographic  conditions 
existing  in  the  region  of  the  Santa  Barbara  Channel  (SBC) 
for  September  2008  are  described  next  to  provide  the  broad 
context  of  the  RaDyO  observations  made  in  the  eastern 
portion  of  the  SBC  from  3-25  September  2008  in  the 
vicinity  of  34.2053°  N,  119.6288°.  Please  note  that 
throughout  the  paper  times  are  given  as  UTC  and  use  the 
year  day  definition  (i.e.,  1  January,  noon  =  year  day  1.5). 
More  detailed  papers  in  this  section  and  elsewhere  present 
more  specific  results  of  the  experiment. 

3.1.1.  Wind  Stress  Patterns  and  Variability 

[78]  Many  different  factors,  which  vary  on  a  broad  range 
of  spatiotemporal  scales,  must  be  considered  to  understand 
and  quantify  the  SBC’s  meteorology,  which  is  complex  and 
thus  demanding  of  numerical  model  simulations  [see 
Carvalho  et  al.,  2012],  The  purpose  of  this  section  is  to 
summarize  the  wind  patterns  and  their  variability  as  the 
wind-forcing  conditions  for  the  SBC  RaDyO  experiment. 
First,  climatological  wind  stress  analyses  for  coastal  Cali¬ 
fornia  were  examined  and  satellite  wind  stress  data 
(QuikSCAT)  were  used  to  document  the  wind  fields  over  the 
waters  offshore  of  central  and  southern  California  (Figure  5). 
Time  series  meteorological  and  atmospheric  flux  records 
were  obtained  from  the  R/V  Kilo  Moana,  R/P  FLIP,  and  the 
NDBC  buoy. 


[79]  The  wind  pattern  based  on  QuikSCAT  satellite  wind 
stress  data  from  5  to  28  September  were  quite  regular  with 
winds  persistently  directed  southeastward  along  the  Cali¬ 
fornia  coast  and  thus  almost  always  upwelling  favorable 
(Figure  5).  This  wind  pattern,  caused  by  relatively  high 
atmospheric  pressure  (anticyclone)  in  the  eastern  North 
Pacific  and  relatively  low  (thermally  induced)  pressure  over 
the  southwest  U.S.  in  the  late  summer  and  early  autumn 
[e.g.,  Dorman  and  Winant,  1995;  Winant  and  Dorman, 
1997],  varies  on  the  synoptic  scale  with  strengthening  and 
weakening  of  the  pressure  systems  and  shifts  in  their  posi¬ 
tion.  No  major  low-pressure  systems  passed  through  the  site 
during  the  experimental  period,  which  is  consistent  with 
seasonal  climatology  [Winant  and  Dorman,  1997]. 

[so]  Unfortunately,  no  satellite  wind  data  are  available 
within  the  SBC  because  of  nearness  to  land  and  limited 
spatial  resolution.  However,  wind  data  collected  in  the  SBC 
from  R/V  Kilo  Moana,  R/P  FLIP,  and  the  NDBC  buoy  show 
that  the  winds  are  generally  directed  toward  the  east  at  the 
study  site  (Figure  6)  [ Zappa  et  al,  2012;  Carvalho  et  al., 
2012.].  The  meteorological  data  collected  from  R/V  Kilo 
Moana  and  R/P  FLIP  were  virtually  identical  when  the 
platforms  were  near  each  other.  The  observed  change  in 
direction  of  the  wind  from  generally  southeastward  to  east¬ 
ward  in  the  vicinity  of  Point  Conception  is  well  documented 
and  winds  roughly  follow  the  coastline  [see  Dorman  et  al., 
1999;  Skyllingstad  et  al.,  2001],  While  there  is  a  strong 
diurnal  signal  in  the  westerly  wind  speed  in  the  SBC,  only 
very  brief  wind  reversals  were  observed  at  the  RaDyO  site 
during  the  period  of  21-23  September  (YD  265-267)  as 
indicated  in  Figure  6.  These  reversals  may  be  indicative  of 
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Figure  5.  Satellite  wind  vectors  (from  QuikScat)  and  sea  surface  temperature  (from  MODIS,  daytime) 
images  for  the  periods  of  (a)  6-12  September  2008  (YD  249-256),  (b)  13-20  September  2008  (YD 
257-264),  and  (c)  21-28  September  2008  (YD  265-272). 
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Figure  6.  (a)  Wind  speed  time  series  data  obtained  from  the  R/V  Kilo  Moana  (blue)  and  R/P  FLIP  (red). 
Mean  and  standard  deviations  were  computed  using  R/V  Kilo  Moana  data,  (b)  Similar  figure  for  wind 
direction  with  same  color  coding. 


atmospheric  eddies  which  sometimes  set  up  over  the  region 
or  may  possibly  be  associated  with  atmospheric  Kelvin 
waves  [e.g.,  Kessler  and  Douglas,  1991;  Dorman  and 
Winant,  1995,  2000].  The  atmospheric  modeling  work  of 
Carvalho  et  al.  [2012]  provides  more  insights  into  the 
atmospheric  dynamics. 

[si]  As  indicated  by  time  series  of  wind  speed  data  shown 
in  Figure  6a,  three  fairly  well-defined  periods  of  mid-to- 
high-to-low  wind  speeds  occurred  during  the  RaDyO  SBC 
experiment  with  average  wind  speeds  (1)  for  9-15  Septem¬ 
ber  (YD  253-259)  of  4.9  m/sec,  (2)  for  16-19  September 
(YD  260-263)  of  7.6  m/sec,  and  (3)  for  20-23  September 
(YD  264-267)  of  3.7  m/sec  (Figure  6a).  These  three  periods, 
which  are  associated  with  the  relative  strengthening  and 
weakening  of  the  larger-scale  atmospheric  pressure  systems 
over  the  Pacific  and  southwest  U.S.,  motivate  further 
inspection  of  the  oceanographic  responses.  The  greatest 
wind  speeds  (~12  m/sec)  for  the  study  occurred  during  the 
afternoon  of  16  September  (YD  260).  Diurnal  variations  in 
wind  speed  are  clearly  evident  in  the  time  series  shown  in 
Figures  6a  with  peak  winds  typically  occurring  in  the  early 
afternoons.  Figure  7  shows  the  summary  of  the  atmospheric 


conditions  for  the  Santa  Barbara  Channel  experiment.  The 
strong  diurnal  cycle  in  the  wind  speeds  and  wave  height 
throughout  the  experiment  is  confirmed  in  the  spectra  of  the 
wind  speed,  Ui0,  and  significant  wave  height,  FIS,  shown  in 
Figure  8.  The  mean  wind  speed  for  the  experiment  was 
5.2  m/sec  with  a  standard  deviation  of  2.8  m/sec,  typical  for 
this  time  of  year  [e.g.,  Winant  and  Dorman ,  1997;  Dorman 
and  Winant,  2000].  The  data  reported  here  have  been  used 
by  Carvalho  et  al.  [2012]  for  comparison  with  their  atmo¬ 
spheric  model  simulations,  which  elucidate  the  dominant 
regional  atmospheric  processes  for  the  RaDyO  experiment. 
3.1.2.  Oceanographic  Conditions 
[82]  The  Santa  Barbara  Channel  RaDyO  field  experiment 
focused  on  a  relatively  small  geographical  region  since  the 
scales  of  interest  were  primarily  on  order  of  a  few  hundred 
meters  in  the  horizontal  scale  and  up  to  a  few  weeks  in 
duration.  Flowever,  it  was  important  to  characterize  larger- 
scale  phenomena  that  affect  the  smaller-scale  processes  as 
well.  In  addition,  the  more  localized  RaDyO  observations 
have  been  used  to  test  atmospheric  and  oceanographic 
models  devoted  to  the  larger  scales  [Carvalho  et  al.,  2012; 
Dong  et  al.,  2012],  With  this  motivation,  we  provide  a  brief 
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YD  — * —  Short-Wave  Radiometer  -  up  (outgoing) 

— • —  Short-Wave  Radiometer  -  down  (incoming) 
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Figure  7.  Environmental  conditions  during  the  Santa  Barbara  Channel  Experiment,  September  2008, 
recorded  from  R/P  FLIP  (Melville  et  al.,  manuscript  in  preparation,  2012).  Note  the  strong  diurnal  cycle 
in  the  wind  speed  and  wave  height  while  the  wind  direction  remains  relatively  fixed  with  only  a  few 
episodes  of  changing  direction. 


17  of  39 


C00H20 


DICKEY  ET  AL.:  INTRODUCTION 


C00H20 


Figure  8.  Spectra  of  30  min  averaged  wind  speed  Uio  and  significant  wave  height  Hs  for  the  entire  Santa 
Barbara  Channel  experiment  record,  in  September  2008,  showing  a  clear  diurnal  cycle  (Melville  et  ah, 
manuscript  in  preparation,  2012). 


summary  of  the  oceanographic  conditions  for  the  period  of 
the  SBC  RaDyO  experiment  in  the  context  of  the  three 
meteorological  forcing  regimes. 

3. 1.2.1.  Currents 

[83]  Horizontal  current  data  for  the  SBC  experiment  were 
derived  from  three  sources:  high-frequency  (HF)  radar  using 
stations  located  along  the  coast  for  near-surface,  channel¬ 
wide  currents,  ADCP  data  collected  from  R/P  FLIP,  R/V 
Kilo  Moana,  and  two  AUVs  (Figure  4).  The  high-frequency 
radar  surface  current  data  (Figure  9)  characterize  the  evolv¬ 
ing  spatial  patterns  of  near-surface  currents  (averaged  over 
~lm  depth),  which  are  important  for  determining  advection 
of  substances  in  surfactants  and  organic  and  inorganic  par¬ 
ticulate  and  dissolved  matter,  all  of  which  affect  optical 
properties.  In  addition,  the  surface  wavefields  are  affected  by 
the  regional  surface  currents  and  vice  versa.  The  uplooking 
ADCP  data  collected  from  R/P  FLIP  (~  1  —25  m)  provide 
high-resolution  current  data  at  depths  of  primary  interest  for 
the  experiment.  Finally,  the  R/V  Kilo  Moana  ADCP  data 
(~  10-200  m)  complete  the  current  records  for  the  upper 
ocean.  There  are  a  few  time  gaps  in  the  R/P  FLIP  data  as 
seen  in  Figure  10.  Evening  gaps  in  the  reported  R/V  Kilo 
Moana  ADCP  data  arose  from  necessary  excursions  of  more 
than  1.5  km  from  the  R/P  FLIP  mooring  site.  These  are  the 
most  problematic  because  of  severe  aliasing  of  current  data 
accentuated  primarily  by  the  diurnal  wind-forcing  and  sec¬ 
ondarily  the  tides.  For  these  reasons,  the  subsurface  current 
records  need  to  be  interpreted  with  these  deficiencies  in 
mind. 

[84]  We  begin  the  discussion  with  the  HF  radar  current 
maps  shown  in  Figure  9,  which  are  a  subset  (12,  19,  and 
23  September;  YD  256,  263,  and  267)  of  the  entire  record. 
Data  shown  in  these  figures  are  centered  on  the  time  periods 
of  the  three  wind-forcing  regimes  (Figure  6)  of  daily  HF 
radar  current  maps,  which  are  available  on  the  Website 
http://www.icess.ucsb.edu/iog/archive/25_lTi-_means_menu.php. 


Currents  using  HF  data  were  computed  using  25  h  averages 
to  minimize  tidal  effects.  It  is  worth  noting  that  currents 
measured  at  depths  of  even  a  few  meters  by  ADCP  current 
profilers  can  be  quite  different  from  HF -measured  (~ 1  m) 
surface  currents  [e.g.,  Siegel  and  Lohrmann,  2010],  These 
differences  are  not  unexpected  considering  near-surface 
veering  due  to  Ekman  veering  effects  as  well  as  differences 
in  methodologies,  fundamental  sampling  and  averaging 
schemes,  and  the  lack  of  perfect  colocation  of  all  instru¬ 
ments.  Overlapping  ADCP  current  data  collected  by  the  R/P 
FLIP  and  the  R/V  Kilo  Moana  were  generally  consistent 
when  the  R/V  Kilo  Moana  was  stationed  within  about 
1 .5  km  of  the  moored  R/P  FLIP.  Throughout  the  experiment, 
near-surface  and  upper  ocean  currents  have  a  large  diurnal 
component  associated  with  the  diurnal  wind-forcing.  The 
mean  currents  for  each  of  the  three  designated  periods  are 
usually  several  times  smaller  than  peak  currents  associated 
with  the  diurnal  wind-forcing.  The  flow  pattern  for  the 
period  of  9-15  September  (YD  253-259)  is  characterized  by 
a  cyclonic  eddy  centered  near  the  deepest  portion  of  the  SBC 
(Figure  9a  showing  YD  256  or  12  September).  The  eddy 
extended  nearly  across  the  Channel  (over  20  km  north-south 
and  east-west).  Eastward  flow  was  seen  north  of  Santa  Rosa 
and  Santa  Cruz  Islands.  The  mean  surface  currents  at  the 
RaDyO  site  for  the  first  period  were  3.3  cm/sec  eastward 
and  4.1  cm/sec  southward.  Currents  in  the  upper  20  m 
indicate  some  shifting  of  directions  as  shown  in  Figure  10 
(note  transition  from  generally  northward  to  southward 
around  YD  257).  The  currents  of  the  upper  20  m  were  fairly 
coherent  with  depth,  although  some  vertical  shear  of  hori¬ 
zontal  currents  was  evident. 

[85]  The  greatest  wind-forcing  for  the  experiment 
occurred  during  the  next  period  (16-19  September;  YD 
260-263).  The  eddy  over  the  proximate  center  of  the  basin 
(somewhat  to  the  west  of  the  geometric  center  at  34.25  N, 
120. 1W)  remained  during  the  period;  however,  the  flow  to 
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Figure  9.  HF  surface  current  data  sequence.  Only  data  collected  on  (a)  12  September  2008  (YD  256), 
(b)  19  September  2008  (YD  262),  and  (c)  23  September  2008  (YD  266)  are  shown  here.  Complete  daily 
sequence  may  be  viewed  on  website  www.opl.ucsb.edu/radyo.  The  green  symbol  marks  the  location  of 
R/P  FLIP  during  the  experiment  and  the  yellow  triangles  indicate  the  locations  of  PIF  radar  sites. 
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Figure  10.  Contour  plots  of  currents:  (top)  U  and  (bottom)  V,  both  in  cm/sec  for  depths  of  ~ 1-25  m 
using  the  SIO  uplooking  ADCP. 


the  north  of  Santa  Rosa  and  Santa  Cruz  Islands  intensified 
(Figure  9b).  The  mean  surface  currents  were  9.3  cm  east¬ 
ward  and  2.2  cm/sec  northward.  The  increased  mean  speeds 
and  directional  shifts  may  have  been  caused  by  the  signifi¬ 
cant  increase  of  wind  speeds.  Unfortunately,  there  is  a  data 
gap  in  the  R/P  FLIP  ADCP  record  on  18  September  (YD 
262).  However,  currents  in  the  upper  20  m  appear  to  be 
fairly  coherent  with  depth. 

[86]  The  final  period  of  the  experiment,  20-23  September 
(YD  264-267)  was  marked  by  the  weakest  wind-forcing. 
The  surface  current  patterns  suggested  a  transitional  flow 
regime  as  the  flow  changed  direction  from  generally  east¬ 
ward  to  westward  through  the  east  end  of  the  SBC 
(Figure  9c).  The  flow  pattern  in  the  vicinity  of  the  experi¬ 
mental  site  on  22  September  (YD  266)  shows  very  weak 
surface  currents;  however  fairly  strong  northward  surface 
currents  developed  to  the  east  of  the  site.  Inspection  of  the 
complete  sequence  of  current  images  suggests  the  possibility 
of  a  submesocale  eddy  developing  and  propagating  west¬ 
ward  in  the  eastern  portion  of  the  SBC.  The  mean  surface 
currents  at  the  site  for  the  final  period  were  2.5  cm/sec 
eastward  and  4.4  cm/sec  northward,  consistent  with  the  flow 
pattern  of  Figure  9c.  Perhaps  because  of  the  transitory  sur¬ 
face  currents,  the  directions  of  currents  at  the  RaDyO  site 
were  highly  variable:  southward  on  20  September  (YD  264), 
near  zero  on  21  September  (YD  265),  strong  toward  the 
north  on  22  September  (YD  266),  and  strong  toward  the 
north-northeast  on  23  September  (YD  267).  Again,  the  R/P 
FLIP  ADCP  data  indicate  coherent  current  structure  with 
depth  (Figure  10).  R/P  FLIP  and  other  ADCP  data  agree  at 
17  m. 

[87]  As  explained  by  Harms  and  Winant  [1998],  the 
observed  surface  currents  of  the  SBC  are  thought  to  result 
from  the  additive  effects  of  a  larger-than-SBC  scale  flow  and 
a  cyclonic  circulation  within  the  SBC  itself.  They  presented 
synoptic  views  of  SBC  circulation  to  include  the  following 
regimes:  upwelling,  relaxation,  cyclonic,  propagating 
cyclones,  flood  east,  and  flood  west.  The  HF  radar  data 
collected  for  our  experiment  appear  to  best  conform  to  the 
cyclonic  model  and  perhaps  the  propagating  cyclones  model 
toward  the  end  of  the  experiment.  Finally,  the  strong  diurnal 
(land-sea  breeze)  forcing  and  tides  played  significant  roles  in 


the  current  variability  within  the  SBC  and  at  the  site  as  well. 
More  details  concerning  the  circulation  within  the  SBC  are 
presented  by  Dong  et  al.  [2012]. 

3. 1.2.2.  Distributions  and  Evolutions  of  Physical  and 
Optical  Properties 

[88]  The  distributions  of  physical  and  optical  properties 
and  their  evolution  during  the  SBC  experiment  are  consid¬ 
ered  next  by  examining  available  satellite  images  and  the 
R/V  Kilo  Moana  and  R/P  FLIP  data  collected  at  the  SBC 
RaDyO  site.  One  of  the  challenges  is  to  distinguish  the  var¬ 
iability  at  the  site  according  to  localized  forcing  versus 
advection.  Thus,  meteorological  and  current  data  are  needed 
to  assist  in  interpretation.  Clearly,  fully  3-dimensional  mod¬ 
eling  of  the  SBC  for  the  experiment  will  enable  more  detailed 
testing  of  hypotheses  and  identifying  the  most  relevant 
processes  and  quantifying  their  scales  [e.g.,  Dong  et  al., 
2012]. 

[89]  As  discussed  earlier,  the  large-scale  forcing  and 
oceanographic  response  along  the  northern  and  central 
California  coast  during  our  experiment  was  dominated  to  a 
large  degree  by  generally  along-coast  (generally  toward  the 
southeast)  upwelling  favorable  winds.  This  is  evident  in  the 
MODIS  sea  surface  temperature  and  QuikSCAT  wind  sat¬ 
ellite  data  shown  in  Figure  5.  Wind  shadowing  by  landmass 
extends  southeastward  of  Point  Conception  and  explains 
much  of  the  variability  in  the  sea  surface  temperature  maps 
with  warmer  waters  extending  along  the  California  coast 
south  and  east  of  the  SBC.  In  fact,  the  SBC  is  a  transitional 
region  where  cool  upwelled  waters  often  enter  the  west  end 
of  the  SBC  and  warm  waters  move  into  it  from  the  east  end. 
Consequently,  the  SBC  is  typically  characterized  by  com¬ 
plicated  flows  and  hydrography  as  well  as  optical  property 
distributions  with  the  occurrence  of  eddies  and  fronts  of 
varying  scales  as  discussed  in  the  previous  section. 

[90]  The  site  of  RaDyO,  which  is  in  the  eastern  portion  of 
the  SBC,  appears  to  have  been  in  a  relatively  quiescent 
location  with  respect  to  currents.  However,  about  one  third 
of  the  way  through  the  experiment,  local  peak  daily  winds 
roughly  doubled.  Examination  of  the  HF  radar  surface  cur¬ 
rent  maps  (Figure  9)  also  suggests  a  change  in  the  flow 
regime  at  the  site  with  apparent  intensification  of  a  current 
north  of  the  Channel  Islands.  Again,  this  flow  appears  to 
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have  extended  to  the  RaDyO  site  and  may  partially  explain 
some  of  the  observed  changes  in  physical  and  optical  prop¬ 
erties  described  below. 

[91]  Time-depth  contours  of  temperature,  salinity,  density 
(crt),  and  stratification  are  shown  in  Figure  11.  The  upper 
layer  temperatures  were  relatively  warm  and  the  mixed  layer 
depth  was  around  15  m  through  about  17  September  (YD 
260).  There  was  significant  diurnal  variation  in  the  mixed 
layer  depth  (MLD)  as  forced  primarily  by  a  strong  land-sea 
breeze  signal  and  secondarily  the  daily  cycle  of  surface 
heating  and  cooling.  The  thermocline  extended  downward 
from  the  base  of  the  mixed  layer  by  about  15-20  m  for  the 
first  period  of  wind-forcing.  The  strong  winds  initiated 
around  16  September  (YD  260)  appear  to  have  resulted  in 
deepening  of  the  mixed  layer  to  about  25  m,  cooling  of  the 
surface  layer  waters,  and  sharpening  of  the  thermocline  and 
pycnocline  for  about  2  days  after  which  the  mixed  layer 
shoaled  to  values  of  —  10  m,  even  less  than  those  observed  at 
the  beginning  of  the  experiment. 

[92]  The  general  optical  properties  of  the  SBC  have  been 
reviewed  and  analyzed  in  considerable  detail  by  Toole  and 
Siegel  [2001]  and  Kostadinov  et  al.  [2007].  Optical  proper¬ 
ties  of  the  SBC,  which  are  affected  by  chlorophyll  a,  are 
correlated  to  a  large  extent  with  eastward  advection  of 
upwelled  (cooler  (Figure  5),  nutrient  rich)  waters  from  off 
Point  Conception  as  suggested  by  the  satellite  color  images 
(MODIS)  shown  in  Figure  12.  South  and  east  of  the 
Southern  California  Bight  (SBC)  and  generally  within  the 
more  wind- sheltered  portion  of  the  SCB,  much  lower  values 
of  chlorophyll  a  and  warmer  waters  prevail,  so  when  flow 
through  the  SBC  is  influenced  more  by  westward  flow,  the 
chlorophyll  a  levels  decrease.  The  SBC  is  in  essence  a  bio¬ 
logical  as  well  as  physical  transition  region  between  cooler, 
more  nutrient  rich,  and  more  biologically  productive  waters 
coming  from  the  Point  Conception  area  to  the  west  and  the 
warmer,  less  productive  waters  of  the  SCB  to  the  south  and 
east.  There  are  often  east-west  and  north-south  gradients  in 
chlorophyll  a  and  temperature  because  of  the  prevailing 
current  patterns.  Although  the  spatial  resolution  of  the  sat¬ 
ellite  color  images  shown  here  is  not  great,  there  do  appear 
to  be  larger  concentrations  of  surface  chlorophyll  a  within 
the  eddy  centered  near  the  deepest  portion  of  the  SBC  and  to 
the  west  of  the  RaDyO  study  site  (Figure  12).  Similar 
observations  of  higher  chlorophyll  a  levels  associated  with 
such  basin-scale  eddies  that  can  extend  to  at  least  200  m 
depth  have  been  reported  by  Nishimoto  and  Washburn 
[2002],  Anderson  et  al  [2006],  and  Anderson  et  al.  [2008]. 
It  is  thought  that  these  eddies  are  convergent  and  lead  to  high 
levels  of  chlorophyll  a  and  fish  stocks.  There  appear  to  be 
submesoscale  eddies  within  the  SBC  as  well  [Beckenbach 
and  Washburn,  2004;  Dong  et  al,  2012]  and  these  too 
may  have  a  major  influence  on  chlorophyll  a  distributions 
and  optical  properties  in  general  [e.g.,  Bassin  et  al.,  2005]. 

[93]  Time  series  of  dissolved  oxygen,  chlorophyll  a,  and 
beam  attenuation  coefficient  are  displayed  as  functions  of 
depth  in  Figure  11.  Dissolved  oxygen  is  relatively  uniform 
in  the  upper  40-65  m  with  some  apparent  influence  by  the 
increased  winds  during  the  period  of  16-20  September  (YD 
260-264).  The  maximum  vertical  gradient  in  dissolved 
oxygen  lags  that  of  the  MLD  by  about  1  day.  The  chloro¬ 
phyll  a  maximum  layer  generally  resides  just  below  the 
mixed  layer,  but  there  are  some  significant  departures.  The 


most  impressive  feature  occurs  between  depths  of  about  10 
and  40  m  around  20-22  September  (YD  264-266).  The 
origin  of  this  high  chlorophyll  a  feature  is  not  certain.  It  is 
possible  that  it  may  have  resulted  from  advection  as  there 
was  a  change  in  current  patterns  around  this  time.  However, 
another  possibility  is  that  the  prior  wind  event  may  have 
been  sufficient  to  cause  entrainment  of  some  of  the  subsur¬ 
face  chlorophyll  a  maximum  waters  and  more  nutrient  rich 
waters  from  depth,  thus  possibly  stimulating  a  phytoplank¬ 
ton  bloom  with  a  time  lag  of  a  couple  of  days. 

[94]  The  beam  attenuation  coefficient  data  show  higher 
values  within  and  just  below  the  mixed  layer,  but  these  fall 
off  quickly  below  for  the  first  two  periods  of  the  experiment. 
However,  similar  changes  as  those  noted  for  chlorophyll  a 
occur  around  20-22  September  (YD  264-266).  In  this  case, 
the  beam  attenuation  coefficient  maximum  lies  very  near  the 
base  of  the  MLD  whereas  the  chlorophyll  a  maximum  was 
deeper  by  about  10  m.  Whether  this  feature  was  caused 
locally  or  resulted  from  advection  will  be  explored  by  ana¬ 
lyzing  data  collected  by  the  Odyssey’  AUV  and  the  REMUS 
AUV  [Moline  et  al.,  2012]  and  by  using  general  circulation 
models  [e.g.,  Dong  et  al,  2012], 

4.  Background  for  the  Hawaiian  RaDyO  Study 

[95]  The  Hawaiian  RaDyO  study  region  lies  in  the  sub¬ 
tropical  North  Pacific  gyre.  The  clockwise  wind  pattern  of 
the  subtropical  North  Pacific  gyre  has  been  documented  in 
several  studies  [e.g.,  Hellerman  and  Rosenstein,  1983].  The 
specific  region  of  the  Hawaiian  RaDyO  study  was  approxi¬ 
mately  17.5  to  18.0  N,  155.5  to  160.0  W  (Figure  13).  The 
wind  patterns  here  generally  vary  from  northeasterly  to 
easterly  (Figure  14)  and  display  some  seasonality  with 
maximum  trade  winds  occurring  in  summer  months  and 
more  moderate  winds  during  winter. 

[96]  The  clockwise  circulating  oceanic  gyre  is  driven  by 
these  winds.  The  wind  stress  curl  is  particularly  important 
for  the  ocean  circulation  as  demonstrated  by  Munk  [1950]. 
The  wind-driven  currents  in  the  RaDyO  study  region  lie  on 
the  southern  side  of  this  gyre  and  are  part  of  the  North 
Equatorial  Current  (NEC).  The  surface  currents  of  the  NEC 
are  directed  predominantly  from  east  to  west  in  latitudes  of 
roughly  10  to  20  N  and  thus  through  the  RaDyO  study  site 
(-17.5-18  N). 

[97]  The  starting  point  for  the  RaDyO  time  series  obser¬ 
vations  from  R/P  FLIP  and  the  accompanying  R/V  Kilo 
Moana  was  almost  due  south  of  the  southernmost  tip  of  the 
Big  Island.  R/P  FLIP  was  allowed  to  drift  for  over  500  km 
(generally  westward)  in  response  to  both  winds  and  currents 
to  avoid  difficulties  in  mooring  it  in  deep  waters  (Figures  13 
and  14).  The  water  depth  was  nominally  4500  to  5000  m. 
The  choice  of  this  study  region  was  based  upon  desires  (1)  to 
do  observations  in  open  ocean  conditions  with  strong,  per¬ 
sistent  winds  and  long-fetch  surface  gravity  waves,  (2)  for 
sampling  in  clear  (Case  1)  waters,  and  (3)  for  relatively  short 
transit  times  to  the  study  region  for  research  vessels  origi¬ 
nating  from  a  major  port  (i.e.,  Honolulu).  The  R/V  Kilo 
Moana  generally  operated  within  about  1.5  km  of  R/P  FLIP 
once  the  R/P  FLIP  reached  the  experimental  starting  location 
shown  in  Figure  13;  periodic  excursions  away  from  R/P 
FLIP  (few  kilometers)  were  required  for  disposal  purposes. 
The  R/V  Kilo  Moana  began  sampling  at  R/P  FLIP’S 
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Figure  11.  (left)  Time-depth  contour  plots  of  temperature,  salinity,  crt,  and  stratification  using  data  collected  with  the  R/V 
Kilo  Moana  and  MASCOT  CTDs  to  100  m  depth.  Depth  of  the  mixed  layer  based  on  a  0.5°C  criterion  is  also  displayed  in 
white  curves,  (right)  Same  as  Figure  11,  left,  for  dissolved  oxygen  in  ml/1,  chlorophyll  a  in  mg/m\  and  beam  cpg  (at  650  nm) 
in  m_1  to  100  m  depth. 
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Figure  12.  Satellite  color  images  indicating  surface  distributions  of  chlorophyll  a  from  north  of  Point 
Conception  and  within  the  Southern  California  Bight  for  the  periods  of  (a)  5-12  September  2008  (YD 
249-256),  (b)  13-21  September  (YD  257-264),  and  (c)  22-29  September  (YD  265-272). 
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Figure  13.  Map  showing  positions  of  the  R/P  FLIP  and  R/V  Kilo  Moana  during  the  Hawaii  RaDyO 
experiment  from  30  August  to  15  September  2009  (YD  242-258).  R/P  FLIP  was  allowed  to  drift  and 
the  R/V  Kilo  Moana  was  positioned  near  R/P  FLIP  for  most  of  the  experiment. 


approximate  starting  point  (where  it  was  flipped  into  ver¬ 
tical  position)  on  YD  242  (30  August  2009)  and  its  sam¬ 
pling  ended  on  YD  258  (15  September  2009).  R/V  Kilo 
Moana  measurements  made  before  YD  242  (30  August 
2009)  were  done  in  transit  or  during  a  diversion  for  mea¬ 
surements  at  a  location  in  the  lee  of  the  Big  Island.  Primary 
R/P  FLIP  time  series  span  the  period  of  approximately  YD 
244  (1  September  2009)  through  YD  257  (14  September 
2009).  The  REMUS  AUV  was  also  deployed  from  the  R/V 
Kilo  Moana  on  YD  243,  247,  and  248  (31  August  and 
4  and  5  September  2009)  in  order  to  obtain  horizontal 
spatial  variability  data  as  well  as  to  obtain  a  unique  view  of 
the  ocean  environment  [. Moline  et  al,  2012].  Satellite  wind, 
sea  surface  temperature  (Figure  15)  and  ocean  color  data 
(Figure  20)  provided  regional  atmospheric  and  oceano¬ 
graphic  context  for  the  study. 

4.1.  Results  for  the  Hawaiian  RaDyO  Study 

4.1.1.  Atmospheric  Conditions 

[98]  The  atmospheric  conditions  for  the  Hawaii  experi¬ 
ment  were  documented  using  satellite  and  in  situ  measure¬ 
ments  from  the  R/V  Kilo  Moana  and  the  R/P  FLIP  in  a 
similar  fashion  to  that  of  the  SBC  experiment  and  with  the 
same  general  motivations  and  objectives.  First,  we  focus 
only  on  the  wind  conditions  for  context  as  details  concerning 
other  meteorological  observations  are  presented  in  Zappa 
et  al.  [2012]  and  Melville  et  al.  (manuscript  in  preparation, 
2012). 

[99]  The  well-documented,  steady  trade  winds  of  the 
region  prevailed  during  the  Hawaii  RaDyO  experiment  as 
evidenced  by  satellite  QuikSCAT  data  for  the  periods  of  YD 
241-248,  249-256,  and  257-264  (29  August  to  21  Sep¬ 
tember  2009;  Figure  15).  Wind  direction  time  series  data 
collected  from  the  R/V  Kilo  Moana  (Figure  16)  and  R/P 
FLIP  (Figure  17)  also  confirm  the  trade  wind  pattern  with 
winds  generally  from  the  east  (shown  in  Figure  16  as  toward 
the  west  with  a  mean  of  261  deg).  Interestingly,  the  direction 
of  drift  of  R/P  FLIP  began  generally  eastward  with  a  modest 
component  to  the  south  (Figures  13  and  14).  About  halfway 


through  the  drifting  portion  of  the  experiment,  it  took  a 
course  still  to  the  east,  but  with  a  slight  northward  compo¬ 
nent  (Figures  13  and  14).  Careful  inspection  of  the  QuikS¬ 
CAT  images  (Figure  15)  shows  that  there  was  some  spatial 


Figure  14.  (top  to  bottom)  Time  series  of  the  R/P  FLIP 
drift  speed  by  components  (U  and  V),  magnitude,  and  direc¬ 
tion  during  the  Hawaii  experiment.  30  min  averages  were 
used. 
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Figure  15.  Satellite  wind  vectors  (QuikScat)  and  sea  surface  temperature  (MODIS)  images  for  the 
periods  of  (a)  28  August  to  5  September  (YD  241-248),  (b)  6-13  September  (YD  249-256),  and 
(c)  14-21  September  (YD  257-264).  The  westward  trajectory  of  R/P  FLIP  during  the  experiment  is 
indicated  by  the  white  line  and  blue  asterisks. 

variability  in  the  wind  field  over  the  drift  region,  which  may 
account  at  least  in  part  for  the  R/P  FLIP’S  slight  change  in 
direction  as  R/P  FLIP’S  drift  is  controlled  by  both  winds  and 
currents.  The  mean  direction  of  R/P  FLIP’S  drift  was  toward 
269  deg,  which  is  very  close  to  the  mean  wind  speed 
direction  toward  261  deg. 


[too]  Time  series  of  wind  speed  collected  by  R/V  Kilo 
Moana  while  in  the  lee  of  the  Big  Island  and  before  its 
rendezvous  with  R/P  FLIP  clearly  showed  the  island  sha¬ 
dowing  effect  as  did  the  QuikSCAT  wind  data.  The  time 
series  of  wind  speed  and  direction  data  collected  from  R/V 
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Figure  16.  (top  to  bottom)  Time  series  of  R/V  Kilo  Moana  wind  direction  (winds  ‘toward’  convention  is 
used)  and  wind  speed,  surface  temperature,  and  surface  salinity.  Means  for  wind  direction  and  speed  are 
shown  for  five  time  periods  and  for  the  entire  experiment. 


Kilo  Moana  and  R/P  FLIP  during  the  concurrent  sampling 
programs  show  that  although  the  direction  is  quite  steady 
(generally  out  of  the  east  and  toward  the  west),  there  is 
considerable  variability  in  wind  speed  with  peak  winds  of 
about  12  m/sec  occurring  during  YD  246  (3  September) 
as  seen  in  Figure  16.  For  convenience,  we  have  computed 
average  wind  speeds  for  five  periods:  YD  243-245  (31 
August  to  2  September):  9.5  m/sec;  YD  245-248  (2-5 
September):  10.9  m/sec;  YD  248-252  (5-9  September): 
9.2  m/sec;  YD  252-256  (9-13  September):  8.5  m/sec,  and 
YD  256-258  (13-15  September):  5.8  m/sec  (Figure  16).  The 
average  wind  speed  for  the  entire  record  was  8.5  m/sec  with 
a  standard  deviation  of  2.6  m/sec.  After  the  first  few  days, 
wind  speed  increased  and  leveled  off  before  decreasing  and 
generally  leveling  off  again.  Interestingly,  the  peak  wind 
speeds  for  the  Hawaiian  experiment  were  very  similar  to 
those  observed  during  the  SBC  experiment  (compare  wind 
speed  data  in  Figure  16  with  those  in  Figure  6).  However, 
the  experimental  period  wind  speed  average  for  the  Hawaii 
experiment  was  3.3  m/sec  greater  than  that  of  the  SBC 
experiment.  Both  display  diurnal  signals,  though  consider¬ 
ably  less  so  for  the  Hawaiian  experiment.  Figure  17  shows 
the  summary  of  the  atmospheric  conditions  for  the  Hawaiian 
experiment.  Diurnal  cycles  in  the  wind  speeds  and  wave 


heights  are  evident  for  the  duration  of  the  experiment.  More 
detailed  meteorological  discussions  may  be  found  in  Zappa 
et  al.  [2012]  and  Melville  et  al.  (manuscript  in  preparation, 
2012). 

4.1.2.  Oceanographic  Conditions 

[101]  The  principal  scales  of  interest  for  the  Hawaiian 
experiment  were  within  approximately  1.5  km  radial  dis¬ 
tance  of  the  R/P  FLIP.  Next,  we  describe  the  general 
oceanographic  conditions  in  the  study  area  from  approxi¬ 
mately  YD  242  through  YD  258  (30  August  to  15  Septem¬ 
ber  2009).  First,  the  current  results  are  presented  and  then 
the  physical  and  optical  data  are  reported. 

4.1.2. 1.  Currents 

[102]  Again,  R/P  FLIP  drifts  in  response  to  wind  and 
current  forcing.  The  relative  contributions  of  the  two  factors 
cannot  be  directly  determined.  However,  for  the  present 
experimental  region,  the  winds  and  currents  historically  are 
generally  directed  toward  the  west.  R/P  FLIP’S  drift  fol¬ 
lowed  a  westward  course  close  to  that  of  the  mean  winds 
(Figures  13,  14,  15,  and  16).  R/P  FLIP’S  drift  speed  by 
components,  magnitude,  and  direction  are  shown  in 
Figure  14.  About  halfway  through  the  experiment,  there  is  a 
shift  from  a  slight  southward  component  to  a  modest 
northward  component  with  an  experimental  mean  northward 
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Figure  17.  Environmental  conditions  during  the  Hawaii  Experiment,  September  2009,  recorded  from 
R/P  FLIP  (Melville  et  al.,  manuscript  in  preparation,  2012).  Note  the  correlation  between  the  wind  speed, 
friction  velocity  from  the  eddy  flux  measurements,  and  the  wave  height  over  the  course  of  the  experiment. 
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Figure  18.  Contour  time  series  of  R/P  FLIP  ADCP  velocity  data  ((top  to  bottom)  U-component, 
V-component,  and  speed)  in  the  upper  12  m.  All  velocities  are  in  cm/sec. 


component  of  less  than  0.2  cm/sec.  The  mean  drift  speed  for 
the  experiment  is  34.3  cm/sec  toward  269  deg..  This  speed  is 
somewhat  greater  than  historical  reports  of  mean  flow, 
probably  because  of  the  added  effect  of  wind  on  FLIP. 

[103]  Subsurface  currents  were  obtained  from  an  uplook- 
ing  ADCP  mounted  at  25  m  on  the  R/P  FLIP  and  a  down¬ 
looking  ADCP  mounted  on  the  hull  of  the  R/V  Kilo  Moana. 
Platform  drift  is  included  in  the  computation  of  the  reported 
currents  for  both.  The  R/P  FLIP  data  span  approximately  the 
upper  12  m  of  the  water  column.  Current  time  series  for 
depths  of  5  and  10  m  are  directed  generally  westward  with 
some  north-south  oscillations  associated  with  tidal  (primar¬ 
ily  semidiurnal)  and  possibly  near  inertial  motions.  Upper 
ocean  current  speeds  generally  increased  by  a  factor  of 
roughly  1.5  to  2.0  during  the  experiment.  The  current 
direction  was  quite  consistently  toward  the  west,  but  with 
some  variability  related  to  tides  as  was  the  case  for  the  R/P 
FLIP  drift.  The  currents  in  the  upper  12  m  were  generally 
uniform  in  speed  and  direction  as  suggested  by  data  shown 
in  Figure  18. 

[104]  R/V  Kilo  Moana’s  300  kHz  ADCP  current  data 
indicated  some  vertical  shear  in  the  horizontal  currents 
appears  near  the  base  of  the  mixed  layer  and  oscillations 
associated  with  tides  (Figure  19).  The  currents  measured 
during  our  experiment  are  consistent  with  those  reported  by 
previously  investigators,  though  slightly  stronger  than  those 
simulated  using  various  numerical  models.  The  most  inter¬ 
esting  aspect  is  the  shift  in  currents  about  midway  through 
the  experiment,  which  may  result  from  a  combination  of 
variations  in  wind-forcing  and  the  movement  of  our  plat¬ 
forms  into  a  different  water  mass  (described  below)  with  a 
differing  current  regime. 


4.I.2.2.  Distributions  and  Evolutions  of  Physical 
and  Optical  Properties 

[105]  The  spatial  variability  of  sea  surface  temperature  of 
the  region  (roughly  17.5-18.0  N,  155.5-160.0  W)  of  the 
Hawaiian  experiment  is  illustrated  in  the  MODIS  satellite 
SST  data  displayed  in  Figure  15.  These  maps  indicate  hori¬ 
zontal  spatial  variability  with  SST’s  roughly  in  the  range  of 
26-27°C  and  with  slightly  warmer  (by  ~1°C)  waters 
generally  lying  to  the  west  of  approximately  157  W  and  to 
the  south  of  the  study  region.  Satellite-based  (MODIS) 
chlorophyll  a  concentration  distributions  are  shown  for 
approximately  the  same  periods  in  Figure  20.  There  was 
considerable  cloud  cover  negating  observations  of  large 
areas,  however  the  influence  of  the  archipelago  is  evident 
with  higher  chlorophyll  a  concentrations  in  the  wind  lee  of 
the  islands  as  expected  based  on  earlier  experiments  includ¬ 
ing  E-FLUX.  These  satellite  data  reflect  only  near-surface 
conditions  and  subsurface  spatial  distributions  likely  differ. 

[106]  Diurnal  oscillations  in  near-surface  temperature  of 
about  0.1-0.2°C  (peak-to-peak)  are  evident  by  inspecting 
time  series  temperature  records  from  R/P  FLIP  (Figures  16 
and  21).  These  records  also  show  relatively  steady  SST’s 
of  about  26.3°C  from  YD  242  (30  August)  until  YD  250 
(7  September),  increasing  SST’s  from  this  point  until  about 
YD  253  (10  September)  when  SST’s  reach  about  27.3°C, 
with  only  modest  change  until  the  end  of  the  experiment  on 
about  YD  258  (15  September).  The  shift  appears  to  be 
caused  by  regional  water  mass  variability,  as  based  on  the 
satellite  SST  map  data,  opposed  to  local  forcing  at  the 
drifting  R/P  FLIP  although  the  greatest  wind-forcing  of 
the  experiment  occurred  prior  to  YD  250  (7  September)  on 
YD  246  (3  September).  Water  mass  change  is  also  supported 
by  salinity  data  obtained  from  the  R/V  Kilo  Moana  as 
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Figure  19.  R/V  Kilo  Moana  ADCP  data  from  17  to  100  m.  (top  to  bottom)  U-component,  V-component, 
speed,  and  shear  in  units  of  cm/sec/m.  The  mixed  layer  depth  is  indicated  with  a  white  line. 


salinity  makes  a  transition  to  lower  values  (from  about 
35.10  psu  to  34.85  psu  at  13  m  depth)  as  temperature 
increases  (Figure  16).  Perhaps  the  most  compelling  evi¬ 
dence  for  the  water  mass  transition  hypothesis  is  the  TS 
diagram  (for  depths  of  0-200  m)  shown  in  Figure  22.  This 
figure  clearly  shows  the  transition  from  cooler  saltier  water 
masses  to  warmer  fresher  water  masses  from  YD  242-250 
(30  August  to  7  September)  to  YD  253-258  (10-15  September). 
Future  model  simulations  should  help  to  sort  out  the  advec- 
tive  versus  local  effects. 

[107]  The  temperature  stack  plot  (Figure  21)  of  time  series 
of  temperature  obtained  from  sensors  mounted  on  R/P  FLIP 
from  the  surface  to  85  m  indicates  increased  stratification 
from  YD  250-253  (7-10  September)  as  well  as  some  vari¬ 
ability  at  depth,  likely  related  to  semidiurnal  internal  tides, 
which  would  be  consistent  with  the  semidiurnal  variability 
noted  in  the  current  records.  Time-depth  contour  plots  of  the 
physical  and  optical  data  are  illustrated  in  Figure  23.  These 
plots  show  some  variability  with  interesting  subsurface 
transitions  occurring  around  YD  250-253  (7-10  September) 
with  a  shoaling  of  the  mixed  layer  followed  by  deepening. 
Changes  in  the  structures  and  magnitudes  of  subsurface 
chlorophyll  a  and  beam  attenuation  coefficient  are  also 


evident.  Since  the  water  mass  differences  are  readily  appar¬ 
ent  before  and  after  the  transition  period  of  YD  250-253  (7— 
10  September),  we  have  also  plotted  in  Figure  24  the  vertical 
structure  of  temperature,  salinity,  density,  chlorophyll  fluo¬ 
rescence,  beam  attenuation  coefficient,  and  dissolved  oxy¬ 
gen  as  average  profiles  for  the  periods  of  YD  242-250 
(30  August  to  7  September)  and  YD  253-258  (10-15  Sep¬ 
tember).  This  figure  provides  a  means  of  summarizing  the 
apparent  water  mass  transition  in  the  upper  200  m  as  fol¬ 
lows:  (1)  increases  in  temperature,  (2)  decreases  in  salinity, 
(3)  decreases  in  density,  (4)  increases  in  stratification,  (5) 
decreases  in  chlorophyll  a  to  about  150  m,  (6)  decreases  in 
beam  attenuation  coefficient  from  about  20  to  80  m  and 
modest  increases  from  about  150-200  m,  and  (7)  a  down¬ 
ward  shift  in  the  depth  of  the  dissolved  oxygen  maximum 
layer. 

5.  Preview  of  Papers  in  the  Special  Section 

[108]  Contributions  to  this  special  section  are  introduced 
categorically  for  convenience  as  follows:  (1)  observational 
papers  based  primarily  on  results  from  the  RaDyO  program, 
(2)  modeling  papers  done  in  conjunction  with  RaDyO,  and 
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Figure  20.  Satellite  color  images  (MODIS)  indicating  sur¬ 
face  distributions  of  chlorophyll  for  the  Hawaiian  study 
region  for  the  periods  of  (a)  28  August  to  5  September 
(YD  241-248),  (b)  6-13  September  (YD  249-256),  and 
(c)  14-21  September  (YD  257-264).  The  westward  trajec¬ 
tory  of  R/P  FLIP  during  the  experiment  is  indicated  by  the 
white  line  and  blue  asterisks. 


(3)  papers  relevant  to  the  theme  of  the  special  section  that 
were  submitted  by  authors  independently  of  RaDyO. 

5.1.  RaDyO  Observational  Papers 

[109]  One  of  the  goals  of  RaDyO  was  to  comprehensively 
observe  and  document  the  environmental  conditions  under 
which  RaDyO  optical  observations  were  made  as  described 
above.  Such  observations  are  necessitated  in  part  by  the 
desire  to  model  radiative  transfer  across  the  air-sea  bound¬ 
ary.  To  do  so  requires  the  detailed  characterization  of  the  sea 
surface  roughness  and  sea  surface  slopes  on  the  microscale. 
To  this  end,  specialized  atmospheric  and  sea  state  observa¬ 
tions  including  air-sea  fluxes  and  surface  topography 
(including  slopes)  for  the  RaDyO  field  experiments  in  the 
Santa  Barbara  Channel  and  in  the  open  ocean  waters  south 


of  Hawaii  are  described  by  Zappa  et  al  [2012]  and  Melville 
et  al.  (manuscript  in  preparation,  2012).  Zappa  et  al.  [2012] 
have  described  and  utilized  new  measurement  capabilities 
for  surface  roughness  and  report  physical  data  pertinent  to 
basic  optical  distortion  processes  related  to  the  air-sea 
interface.  Highlights  of  their  work  include  the  application  of 
a  new  polarimetric  imaging  camera;  this  technique  has 
revealed  a  complex  interaction  of  wind  and  surface  currents 
affecting  surface  roughness.  Additionally,  they  have  repor¬ 
ted  that  breaking  crest  length  spectral  density  is  modulated 
by  the  development  of  the  wavefield  (i.e.,  wave  age)  as  well 
as  the  alignment  of  wind  and  surface  currents  at  the  scale  of 
wave  breaking.  Melville  et  al.  (manuscript  in  preparation, 
2012)  describe  the  impact  of  wave  breaking  and  near-surface 
turbulence  on  imaging  through  the  ocean  surface,  finding 
that  laboratory  simulations  of  these  processes  are  represen¬ 
tative  of  field  data  from  the  three  field  experiments. 

[110]  RaDyO  investigators  report  several  advances  in 
capabilities  to  observe  optical  variability  as  affected  by 
diverse  factors  and  quantified  through  several  variables  on 
fine  space  and  time  scales  as  well  as  spectral  resolution. 
Several  examples  of  these  advances  and  their  application  to 
near-surface  optical  variability  follow. 

[111]  Twardowski  et  al  [2012]  utilized  their  recently 
developed  MASCOT  system  to  measure  the  effects  of  par¬ 
ticles  and  bubbles  on  the  optical  volume  scattering  function 
(VSF)  in  the  surf  zone  off  Scripps  Pier  during  the  initial 
RaDyO  field  exercise.  They  report  how  the  optical  phase 
function  was  affected  by  passing  suspended  sediment 
plumes  and  wave-injected  bubble  plumes.  Measured  VSFs 
were  used  to  determine  particle  size  distributions  and  com¬ 
positions.  Their  inversion  method  results  for  these  data  sets 
showed  qualitative  agreement  with  concurrent  acoustical 
measurements  of  bubbles  and  aggregate  particle  size 
distributions. 

[112]  Unique  measurements  of  bubble  size  distributions 
using  simultaneous  acoustical  and  optical  instrumentation 
are  further  described  by  Czerski  et  al.  [2011].  While  acous¬ 
tical  measurements  of  larger  bubbles  have  been  made  during 
many  studies,  much  less  work  has  been  done  to  measure 
bubbles  with  radii  smaller  than  about  30  pm.  These  are 
thought  to  be  extremely  important  for  optical  properties  and 
radiative  transfer  in  the  upper  ocean,  but  it  is  not  straight¬ 
forward  to  extend  acoustical  measurements  to  these  smaller 
sizes.  To  our  knowledge,  Czerski  et  al.  [2011]  report  the  first 
use  of  simultaneous  acoustical  and  optical  observations  to 
constrain  bubble  coating  parameters.  Their  principal  find¬ 
ings  are  that  organic  coatings  of  the  bubbles  are  relevant  for 
both  acoustical  and  optical  measurements  and  that  the 
inferred  bubble  coatings  (with  a  thickness  of  around  10  nm) 
are  important  for  optical  inversions. 

[113]  With  growing  interest  in  bubble  influences  on  ocean 
optics,  RaDyO  investigators  studied  the  interplay  between 
upper  ocean  turbulence  and  varying  bubble  size  distributions. 
In  particular,  Vagle  et  al.  [2012]  conducted  simultaneous 
measurements  of  upper  ocean  heat  fluxes,  turbulence,  and 
bubble  size  distributions  from  R/P  FLIP  south  of  Hawaii. 
They  found  that  turbulent  dissipation  rates  were  reduced  by 
up  to  a  factor  of  ten  during  times  of  large  downward  heat 
fluxes  (>400  W  m-2).  Vagle  et  al.  [2012]  conclude  that  such 
reductions  can  cause  decreased  concentrations  of  large 
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Figure  21.  Temperature  time  series  for  data  collected  from  R/P  FLIP  sensors  placed  at  18  depths 
between  the  surface  and  85  m. 


bubbles  and  decreases  of  the  effect  of  the  bubble  field  on 
optical  reflectance  by  up  to  a  factor  of  ten. 

[114]  The  sea-surface  microlayer  (SML)  has  been  studied 
intensively  by  Wurl  et  al.  [20 1 1  ]  in  diverse  oceanic  envir¬ 
onments.  They  have  specifically  examined  the  production 
and  fate  of  transparent  exopolymer  particles  (TEP),  which 
have  been  shown  to  have  significant  influences  on  light 
propagation  near  the  air-sea  interface.  Their  work  reinforces 
consensus  on  the  gelatinous  nature  of  the  SML,  which 


influences  microbial  life,  surface  wave  properties,  and  light 
propagation.  Finally,  they  describe  a  conceptual  model  of 
TEP  cycling,  which  suggests  the  importance  of  abiotic  par¬ 
ticle  aggregation  for  TEP  production  in  the  ocean  and  that 
the  SML  plays  a  major  role  in  the  marine  carbon  cycle. 

[115]  High-frequency  measurements  of  the  underwater 
light  field  from  theUnderwater  Porcupine  Radiometer  Sys¬ 
tem  are  reported  by  Darecki  et  al.  [2011].  This  system, 
which  measures  downwelling  light  at  up  to  1  KHz  using  an 
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Figure  22.  Temperature-Salinity  (T-S)  plot  of  R/V  Kilo  Moana  CTD  data  collected  from  the  surface  to 
150  m.  The  color  coding  indicates  blue  for  data  collected  between  YD  242  and  250,  green  for  data 
collected  between  YD  250  and  253,  and  red  for  the  period  of  YD  253  to  258. 


array  of  23  radiometric  sensors,  was  used  to  determine  that 
irradiance  collectors  of  a  few  millimeters  in  diameter  or  less 
are  required  to  adequately  sample  observed  light  flashes 
resulting  from  wave  focusing.  Darecki  et  al.  [2011]  have 
shown  that  some  of  the  brightest  light  flashes,  which  can 
have  durations  of  milliseconds  to  tens  of  milliseconds,  can 
be  an  order  of  magnitude  greater  than  time-averaged  irradi¬ 
ance  values.  They  conclude  that  intensities  of  light  flashes 
decrease  rapidly  as  a  function  of  depth  and  that  intensities  at 
longer  visible  wavelengths  are  greater  than  those  at  shorter 
wavelengths.  Related  work  by  Gernez  et  al.  [201 1]  focused 
on  vertical  variability  of  probability  distributions  of  down¬ 
ward  irradiance  near  the  ocean  surface  under  sunny  skies. 
They  utilized  time  series  measurements  of  fluctuations  in  the 
underwater  downward  irradiance  at  532  nm  (green  light)  to 
quantify  the  nature  of  the  probability  distributions  of 
instantaneous  irradiance.  It  was  shown  that  light  flashes  near 


the  ocean  surface  result  in  highly  skewed  and  heavy  tailed 
probability  distributions  because  of  surface  wave  focusing. 
They  suggest  several  probability  distribution  models  that 
may  be  suitable  for  describing  these  skewed  distributions.  At 
greater  depths  (e.g.,  ~10  m  for  their  observational  situation) 
the  probability  distributions  were  found  to  approach  a  sym¬ 
metric  shape. 

[l  16]  During  RaDyO,  Lewis  et  al.  [2011]  and  Bhandari 
et  al.  [2011]  utilized  their  newly  developed  underwater 
cameras  to  examine  full  spherical  and  polarized  downwelling 
radiance  distributions,  respectively.  The  challenges  facing 
these  two  groups  of  investigators  included  the  need  to  mea¬ 
sure  light  intensify  changes  in  response  to  rapidly  changing 
atmospheric  and  surface  wave  conditions  over  several  orders 
of  magnitude  in  dynamic  range.  The  measurements  made  by 
Lewis  et  al.  [201 1]  are  noteworthy  and  valuable  in  that  their 
fully  resolved  radiance  data  can  be  used  to  derive  all  apparent 
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Figure  23.  (left)  Contour  time-depth  plots  of  temperature,  salinity,  density  (crt),  and  stratification,  (dert/dz)  to  200  m  and 
(right)  dissolved  oxygen,  chlorophyll  a,  and  beam  attenuation  coefficient  at  right. 
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Figure  24.  Average  vertical  profile  data  (0-200  m,  from  R/V  Kilo  Moana  CTD)  for  the  periods  of  YD 
242-250  in  blue  and  YD  253-258  in  red  for  the  following  variables:  (left)  temperature,  salinity,  density 
(sigma-t)  and  (right)  chlorophyll  fluorescence,  beam  attenuation  coefficient,  and  dissolved  oxygen. 


and  inherent  optical  properties  as  well  as  phase  functions.  In 
particular,  they  utilized  their  radiance  data  sets  to  compute 
several  irradiances  (i.e.,  upward  and  downward  vector  and 
scalar  irradiances).  These  quantities  were  further  used  by 
Lewis  et  al.  [2011]  to  determine  other  optically  important 
quantities  such  as  diffuse  attenuation  coefficients.  Impor¬ 
tantly,  they  have  demonstrated  good  agreement  between  their 
data  sets  and  those  obtained  independently  by  other  RaDyO 
investigators.  Bhandari  et  al.  [201 1]  deployed  a  new  camera 
system  called  the  downwelling  polarized  radiance  distribu¬ 
tion  camera  system  (DPOL)  from  R/P  FLIP  in  both  the  Santa 
Barbara  Channel  and  off  Hawaii  during  RaDyO.  The  system, 
which  utilizes  a  combination  of  four  separate  images  simul¬ 
taneously  to  measure  the  Stokes  vector,  measured  the 
downwelling  hemisphere  of  polarized  radiance  distribution 
at  seven  wavelengths  in  the  visible.  This  distribution  is 
valuable  in  that  it  provides  the  most  complete  light  field 
description.  Bhandari  et  al.  [2011]  determined  that  under 
clear  sky  conditions,  the  major  source  of  polarization  is 
refracted  skylight  close  to  the  sea  surface  whereas  deeper  in 
the  water  column  polarization  due  to  light  scattering  both 
increases  with  depth  and  dominates. 

[in]  An  important  theme  of  RaDyO  concerned  under 
water  imaging,  which  was  addressed  by  Chang  and 
Twardowski  [2011].  They  have  computed  a  key  quantity 
relevant  to  underwater  electro-optical  imaging  systems 
called  the  modulation  transfer  function  (MTF).  The  MTF  is 
defined  as  the  Fourier  transform  of  the  optical  point  spread 
function  (PSF).  The  PSF  is  defined  and  discussed  in  detail 
by  Voss  [1991]  who  describes  it  in  physical  terminology  as 


the  image  of  an  unresolved  cosine  source  obtained  from  a 
camera  and  a  point  source  at  opposing  positions.  He  indi¬ 
cates  that  if  were  there  no  scattering  particles  in  the  medium, 
the  PSF  would  be  a  b-function;  however,  with  scattering  the 
image  is  blurred.  Here,  Chang  and  Twardowski  [2011] 
report  that  imaging  parameters  were  correlated  with  the 
characteristics  of  particle  compositions.  They  also  found  that 
optical  properties,  which  were  measured  in  shallow  near 
coastal  waters  (eutrophic)  and  deeper  offshore  (mesotrophic) 
coastal  waters  during  RaDyO,  were  strongly  affected  by 
atmospheric  forcing  conditions. 

[t  is]  Moline  et  al.  [2012]  utilized  autonomous  underwater 
vehicles  (AUVs)  to  map  physical  and  optical  water  vari¬ 
ability  during  RaDyO  in  the  Santa  Barbara  Channel  (SBC) 
and  in  the  open  ocean  off  Hawaii.  Their  measurements  allow 
data  collection  at  scales  and  depths  unattainable  from  con¬ 
ventional  satellite  and  aircraft  remote  sensors.  Data  collected 
from  a  radiometer  mounted  on  an  AUV  during  both  field 
experiments  were  used  to  characterize  optical  constituents 
by  applying  bio-optical  inversion  algorithms.  Their  results 
were  in  good  agreement  with  independent  data  sets  for 
quantities  such  as  optical  absorption. 

5.2.  RaDyO  Modeling  Papers 

[119]  Carvalho  et  al.  [2012]  and  Dong  et  al.  [2012] 
applied  atmospheric  and  oceanic  models,  respectively,  to 
better  define  and  characterize  the  SBC  atmospheric  forcing 
and  oceanic  responses  for  the  period  of  the  RaDyO  experi¬ 
ment.  In  this  section,  Dong  et  al.  [2012]  applied  an  auto¬ 
mated  eddy  detection  scheme  to  a  12  year,  3 -dimensional 
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high-resolution  numerical  oceanic  model  product.  They 
found  a  significant  percentage  of  submesoscale  eddies, 
which  tended  to  be  ageostrophic,  were  generated  around  the 
islands  and  headlands  along  the  coastline.  Three  types  of 
submesoscale  eddies  were  defined  on  the  basis  of  shape.  The 
eddy  fields  characterized  in  this  study  are  valuable  for 
interpreting  optical  data  sets  collected  during  RaDyO  in  the 
SBC. 

[120]  The  next  three  papers  concern  the  challenging 
problem  of  simulating  radiative  transfer  across  the  rapidly 
changing  and  complex  air-sea  interface  and  within  the 
dynamic  and  turbulent  upper  ocean.  Xu  et  al.  [2012]  per¬ 
formed  numerical  simulations  of  radiative  transfer  of  natural 
light  in  turbulent  flows  in  the  upper  ocean.  They  posed  the 
problem  for  turbulent  shear  flow  which  interacted  with  the 
sea  surface  both  with  and  without  surface  waves.  The  phys¬ 
ical  large-eddy  simulations  were  done  for  fluid  motions  and 
transports  of  temperature  and  salinity.  Inherent  optical 
property  (IOP)  variability  was  based  on  turbulent  salinity  and 
temperature  fields.  They  found  that  the  vertical  profile  of  the 
downward  irradiance  is  governed  primarily  by  the  IOP  ver¬ 
tical  structure  with  the  horizontal  variations  being  attributed 
to  the  turbulence.  With  surface  waves,  irradiance  variability 
is  enhanced  via  surface  deformations  caused  by  wave- 
turbulence  interactions.  In  a  related  paper,  Xu  et  al.  [2011] 
report  on  patterns  and  statistics  of  light  polarization  in  the 
water  for  realistic  linear  and  nonlinear  ocean  surface  wave 
conditions.  Direct  simulations  of  linear  and  nonlinear  surface 
waves  using  a  higher-order  spectral  method  set  the  physical 
conditions  for  the  radiative  transfer  model,  which  utilizes  a 
Monte  Carlo  vector  radiative  transfer  simulation.  The  authors 
found  (1)  that  increased  roughness  of  the  surface  leads 
to  smaller  degrees  of  polarization  within  Snell’s  window, 
(2)  variability  of  the  degree  of  polarization  depends  strongly 
upon  surface  roughness,  and  (3)  nonlinearity  of  surface 
waves  primarily  affects  variability  of  in-water  polarizations. 
They  have  validated  their  results  using  RaDyO  data  sets 
collected  during  the  Hawaiian  experiment. 

[121]  The  final  modeling  paper  by  You  et  al.  [2011]  also 
examines  polarized  light  fields  under  the  dynamic  ocean 
surface.  These  authors  developed  a  hybrid  numerical  model 
to  efficiently  compute  and  simulate  polarized  underwater 
light  fields  by  combining  the  advantages  of  3-dimensional 
Monte  Carlo  and  matrix  operator  methods.  This  model  is 
optimized  for  the  dynamic  atmosphere-surface-ocean  system 
in  which  the  temporal  variations  are  limited  to  the  surface 
while  the  atmosphere  and  ocean  can  be  largely  considered  as 
static  within  the  time  frame  of  interest.  This  enables  the 
authors  to  use  a  static  Monte  Carlo  model  for  the  computa¬ 
tionally  intensive  radiative  transfer  calculations  in  the 
atmosphere  and  ocean.  These  calculations  need  to  be  done 
only  once.  On  the  other  hand,  the  interaction  between  the 
light  fields  and  the  ocean  surface  is  highly  dynamic  and  has 
to  be  calculated  for  each  time  instance.  This  interaction  is 
governed  by  the  Fresnel  fonnulae  and  is  rather  straightfor¬ 
ward.  The  static  and  dynamic  parts  of  the  system  are  then 
coupled  using  a  matrix  operator  method.  The  model  uses 
data  sets  (i.e.,  optical  properties  and  surface  wave  slopes) 
collected  during  the  RaDyO  field  experiments.  Model  sim¬ 
ulation  results  for  radiance  and  polarization  fields  were 
found  to  be  qualitatively  consistent  with  their  observed 
counterparts.  The  authors  suggest  that  quantitative 


agreement  would  be  expected  for  colocated  and  synchro¬ 
nized  wave-slope  and  optical  property  measurements.  It  is 
also  implied  that,  given  appropriate  input  data  sets,  model 
simulations  might  be  more  capable  of  providing  a  complete 
set  of  spatial,  temporal  distributions  as  well  as  polarization 
states  of  the  underwater  light  field  comparable  with  most 
advanced  measuring  techniques. 

5.3.  Other  Papers  Relevant  to  the  Theme  of  This 
Special  Section 

[122]  The  final  group  of  papers  were  not  based  on  the 
RaDyO  program,  but  have  relevance  to  the  section  theme. 
The  first  two  by  Salisbury  et  al.  [2011]  and  Gallegos  et  al. 
[2011]  focus  on  regional  studies  of  optical  variability  in 
the  upper  ocean  in  the  outflow  region  of  the  Amazon  River 
(western  tropical  Atlantic)  and  Chesapeake  Bay,  respec¬ 
tively.  The  last  paper  by  Szeto  et  al.  [2011]  considers  optical 
differences  among  four  separate  ocean  basins:  the  Atlantic, 
Pacific,  Indian,  and  Southern  Oceans.  All  of  these  con¬ 
tributions  have  utilized  both  satellite  remote  sensing  and  in 
situ  observations. 

[123]  Salisbury’  et  al.  [2011]  documented  coherences  in 
space  and  time  of  physical  and  optical  properties  both  within 
and  adjacent  to  the  Amazon  River  Plume.  They  report 
freshened,  colored  waters  coming  from  the  Amazon  River 
with  general  coherence  of  salinity  and  colored  dissolved 
matter  absorption  (acdm).  However,  they  note  spatial  pat¬ 
terns  of  these  two  parameters  which  are  not  overlapping. 
Following  the  river  plume’s  trajectory,  they  observed  an 
inverse  relationship  between  salinity  and  acdm,  suggesting 
departure  from  conservative  mixing  for  all  observed  seasons. 
Processes  contributing  to  this  departure,  according  to  the 
authors,  include  areas  of  enhanced  net  primary  productivity; 
they  suggest  the  importance  of  phytoplankton  biomass  or  its 
remineralization  as  contributors  to  the  colored  dissolved 
matter  (cdm)  signal.  They  also  indicate  that  photo-oxidation 
of  cdm  may  be  important  on  longer  time  scales  (days  to 
weeks). 

[124]  Gallegos  et  al.  [2011]  have  utilized  a  long  time 
series  (~25  years)  of  optical  data  collected  in  Chesapeake 
Bay  to  examine  long-term  changes  in  light  scattering  in  that 
body  of  water.  Their  data  sets  come  from  in  situ  Secchi  disk 
and  light  attenuation  records  and  remote  sensing  ocean  color 
satellites.  They  applied  a  bio-optical  model  to  these  data  sets 
to  determine  the  causes  of  observed  trends  in  quantities 
including  the  Secchi  depth,  ZSD,  and  the  diffuse  attenuation 
coefficient  of  photosynthetically  active  radiation  (Kd(PAR)) 
as  well  as  the  product  of  the  two.  Their  model  results  suggest 
a  possible  cause  of  decreasing  backscatter  ratio  to  be 
increasing  proportions  of  organic  detritus  in  the  Bay.  How¬ 
ever,  the  authors  also  indicate  that  increasing  occurrences  of 
large  particle  aggregation  would  be  consistent  with  their 
model  results. 

[125]  The  last  paper  considers  why  the  world’s  oceans 
display  significant  optical  differences.  Szeto  et  al.  [2011] 
have  examined  a  large  volume  of  ocean  color  data  col¬ 
lected  in  situ  and  assembled  by  NASA’s  Ocean  Biology 
Processing  Group  [  Werdell  and  Bailey,  2005].  The  work 
was  stimulated  by  biases  that  appeared  in  the  SeaWiFS 
chlorophyll  algorithm  from  ocean  to  ocean.  The  authors 
were  able  to  rule  out  instrumentation  and  analytical  causes 
(i.e.,  artifacts)  of  the  distinctions  found  among  the  Atlantic, 
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Pacific,  Indian,  and  Southern  Ocean  data  sets.  The  biases  are 
suggested  to  be  attributable  to  regional  differences  in 
absorption  properties  of  detrital  matter  and  phytoplankton. 
Broader  implications  concerning  the  applications  of  ocean 
color  algorithms  to  regions  are  discussed  in  the  authors’ 
conclusions. 
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